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Abstract. Retrievals of a-erosol optical depth (AOD) and related parameters from satel- 
lite measurements typically involve prescribed models of aerosol size and composition, 
and are therefore dependent on how well these models are able to represent the radia- 
tive behaviour of real aerosols. This study uses aerosol volume size distributions retrieved 
from Sun-photometer measurements at 11 Aerosol R.obotic Network (AERONET) island 
sites, spread throughout the world's oceans, as a basis to define such a model for un- 
polluted maritime aerosol Size distributions are observed to be birnodal and approxi- 
mately lognormal although the coarse mode is skewed with a long tail on the low- radius 
end. The relationship of AOD and size distribution parameters to meteorological con- 
ditions is also examined. As wind speed increases, so do coarse-mode volume and radius. 
The AOD and Angstrom exponent (a) show linear relationships with v/ind speed, al- 
though there is considerable scatter in all these relationships, limiting their predictive 
power. Links between aerosol properties and near-surface relative humidity, columnar wa- 
ter vapor, and sea surface temperature are also explored. A reconmiendecl birnodal mar- 
itime model, which is able to reconstruct the AERONET AOD with accuracy of order 
0.01-0.02, is presented for use in aerosol remote sensing applications. This accuracy holds 
at most sites and for wavelengths between 340 nm and 1020 nrn. Calculated lidar ra- 
tios are also provided, and differ significantly from those currently used in Cloud- Aerosol 
Lidar with Orthogonal Polarization (CALIOP) processing. 


1. Introduction 

The size distribution and spectral complex refractive in- 
dex of aerosols are needed to compute properties such as 
their scattering phase function, single scatter albedo, and ex- 
tinction coefficient, which are in turn used to calculate quan- 
tities such as total aerosol optical depth (AOD) from column 
abundance. The information content of measurements from 
current satellite radiometers is insufficient to unambiguously 
retrieve all these parameters, particularly when the (spec- 
tral and directional) behaviour of surface reflectance is un- 
known (Hasekarnp and Landgraf, 2007), For this reason, 
aerosol retrieval algorithms employed by most of these sen- 
sors are required to make assumptions about aerosol micro- 
physical properties and rely on a set of predefined aerosol 
models or components. The assumptions in these aerosol 
retrieval algorithms contribute to differences in retrieved 
AOD, even in the idealised case of a black (non-reflecting) 
surface (Kokhanovsky et a/., 2010). The Polarization and 
Directionality of the Earth’s Reflectance (POLDER) sensor 
is an exception to this, as its measurement capabilities pro- 
vide an increased information content as corap«ared to other 
current sensors [Dubovik ef nl, 2011, Hasekarnp et al, 2011). 

For other sensors, it is therefore of high importance that 
the models used are representative of real aerosol proper- 
ties, The purpose of this study is to develop such a model for 
clean maritime aerosol, using Sun-photometer data from the 
Aerosol Robotic Network (AERONET; Holben et al, 1998). 
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A companion paper, Sayer et at. [2011], describes the appli- 
cation of this model to aerosol retrievals from Sea-viewing 
Wide Field-of-view Sensor (SeaWiFS) measurements. 

Datasets including the optical properties of marine 
aerosol as determined from ground-based measurements, 
aircraft, remote sensing, or theoretical considerations exist, 
and a review of some of these is presented by Smirnov ei al. 
[2002]. In particular, the models Shettle and Fenn [1979] 
(from aircraft measurements) and Gathman [1983] (coastal 
towers, and ships) have been used widely. However, obser- 
vational datasets are typically limited in time and space, 
and differences between the types of instrumentation used 
in these campaigns contribute to significant differences be- 
tween the results {Reid et ai, 2006). Advantages of the 
AERONET data therefore include the opportunity for a 
longer time series, with a wide global distribution, and con- 
sistency between different measurement sites. Such stud- 
ies are also often coastal, such that there may be a non- 
rnaritime component to the aerosol While still a factor 
for AERONET data, this can be minimised through choice 
of remote sites and careful filtering of data. A previous 
AERONET-bfised analysis was performed by Smirnov et al. 
[2003a], although at that time the available data record was 
smaller. 

The aerosol number size distribution d iV(r) /din (r) de- 
scribes the number of aerosol particles with radius in the 
infinitesimal size range r drdln(r). The distribution is also 
sometimes defined as dN(r)/d r, and these two are easily 
related by 

dN{r) _ dln{r)dN(r) ^ IdN’(r) ^ 

dr dr dln{r) r dln(r) ' 

The volume size distribution, calculable for spherical 
aerosol particles as 

dVir) _ 47rr’^ diV(r) 
dln(r) 3 din(r)' 
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describing the aerosol particle volume over the same in- 
finitesinial radius raiige, is also frequently used. The 
AERONET products are defined in terms of the columnar 
volume size distribution and so this convention is adopted 
in the analysis here. The total aerosol columnar particle 
rmmber (Or,) and volume (Cv) are obtained by integrating 
these distributions over all ln(r). 

Frec|uently-used metrics to characteri.se aerosol size dis~ 
tributions include the volume mean radius (?%) as a measure 
of the size of the aerosol particles, where 


Infrv) 



ln(r) 


d V (r) 
dln(?') 


d\n{r) 



dV(r) 

dln(r) 


dln(r) 


(3) 


and the standard deviation (or spread) of the distribution 
(cr) as a measure of the dispersion: 


Note that the spread cr remains the same for both number 
and volume distributions. Hinds [1999] presents some gen- 
eral results for moments of lognormal distributions, includ- 
ing that 
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Section 2 describes the AERONET data used, and prop- 
erties of average size distributions. Section 3 examines the 
effect of meteorology on the size distribution. Next, Sec- 
tion 4 combines the size information with various refrac- 
tive indices to define an average aerosol model which is best 
able to replicate the observed AERONET AODs. Follow- 
ing the definition of this model. Section 5 tests the predic- 
tive powder of relationships observed between wind speed and 
aerosol volume on ship-borne AOD mea.surements, and Sec- 
tion 6 presents calculated Ildar ratios. Finally, Section 7 
summarises the results of the study. 
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In the above the integration is theoretically carried out 
over all ln(r), although in practical applications some min- 
imura and maximum bounds on the radius are defined as 
cutoffs, outside which the aerosol number and volume are 
negligible. The mean radius of the number distribution Vn 
is defined analagously to Equation 3. using dN{r)/dhi{r) in 
place of dV (r) / dln(r) . A third useful quantity is the effec- 
tive radius (ren). the ratio of the third to second moments 
of the number size distribution: 





d N (r) 
d\n{r) 


dln(r) 


d N (r) 
d\n{r) 


dln(r) 


(5) 


The effective radius is more closely related to aerosol ex- 
tinction than the number median radius because scattering 
depends on aerosol cross-sectional area, and distributions 
with similar effective radii (and effective variances, although 
this quantity is not frequently used in aerosol studies) typ- 
ically have similar scattering properties, even if the precise 
modal radii and spreads differ {Hansen and Travis, 1974, 
A'fishchenko et al. 1997). 

Aerosol size distributions are commonly represented as 
a combination of lognormally-distributed components, in 
which case the number size distribution is defined as a sum- 
mation over these (nc) components by 


n i / In(r) - ln(rn,0 

dNjr) ^ ^ Cn,r ^"2 \ cTi 

dln(r) ^ \plr.ai^ 


( 6 ) 


and the modal radius for each component is also its me- 
dian and geometric mean. The equivalent distribution for 
aerosol volume is arrived at by substituting fn with Tv, and 
Cn with Cv. The advantages of iognonnal distributions in- 
clude that their statistical properties are v/ell-known, and 
many available radiative transfer codes are able to take as 
input lognormal distribution paraineters. Fbr individual log- 
normal components, the conversion between the volume and 
number distribution parameters is presented in Appendix A. 


2. AERONET sites and size distribution 
data 

2.1. Sites and data selection criteria 

AERONET data from sites listed in Table 1, and shown 
in Figure 1, are used here to investigate the characteris- 
tics of maritime aerosol. These sites have been chosen due 
to their general remoteness from local sources, to maximise 
the chances of measuring unpolluted maritime aerosol, and 
span a variety of oceans. The stability and pointing ac- 
curacy required to perform the almucantar scans used to 
retrifwe the size distribution means they are impractical 
to perform aboard moving platforms such ships, and so 
island sites represent the closest to open-ocean conditions 
which can be obtained using this technique. Of these sites 
Lanai, Bermuda, and Kaashidhoo were previously studied 
by Smirnov et al. [2003a] . A similar study, also examin- 
ing the effect of wind speed on aerosol properties, w^as per- 
formed for Midway Island by Smirnov et al. [2003b] . The 
ma.in development of this study over previous work is the 
improved data record, in terms of an increased number of 
observations over a larger number of locations, and taking 
advantage of more recent AERONET algorithm improve- 
ments {Holben et al. 2006). Additionally, some meteoro- 
logical aspects, and refractive index, are examined in more 
detail. Ahmad et al. [2010] also used AEiRONET inversions 
to inform aerosol models for satellite retrieval, but with a 
different approach, and did not filter for ‘pure maritime’ 
cases in this way. 

For all sites except Graciosa, only version 2 level 2.0 
(cloud-screened and quality- assured; Smimov et a7, 2000a, 
Holben et ai, 2006) data are ised. Only retrievals from 1999 
onwards are considered as the newer Sun photometers de- 
ployed since then enable a higher data quality. Graciosa is 
a comparatively new AERONET site (established in 2009) 
and some of the earlier measurements are suspected of cloud- 
contamination problems; only level 1.5 (cloud-screened) 
data collected in 2010 are used here as they are thought 
to be largely free of this, although have not yet been for- 
mally raised to level 2.0. Graciosa is used in preference to 
the nearby Azores site (which does have level 2.0 data) as 
the former is in a more remote setting and has a larger data 
record. Despite small data records and/or increased chance 
of influence from other aerosol types, Kaashidhoo, .Ascension 
Island, and Amsterdam Island are included in the analysis as 
the most suitable (or only) sites in their respective regions. 
In particular, Kaashidhoo samples air nuisses transported 
from the Indian subcontinent, south-east Asia, the Arabian 
Sea, and the Southern Indian Ocean {Lobert and Harris, 
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2002), while Ascension Island can be affected by transported 
African biomass burning emissions (Galaniei' et ai, 2000). 
Crozet Island has a small data record, due primarily to fre- 
quent cloud cover, and has the highest elevation above sea 
level (221 m, so still within the marine boundary layer), but 
is included nonetheless as, unlike the majority of other sites, 
it occupies a cool-sea and high- wind environment. 

The AERONET inversion algorithm used to retrieve the 
aerosol size distribution (in 22 logarithmically-spaced size 
bins) and refractive index from Sun-photometer measure- 
ments is described by Dubovik and King [2000] and Dubovik 
et al [2006]. It takes as input diffuse-sky radiances at 440 
nm, 675 nm, 870 nrn, and 1020 nm in the solar almucantar, 
as well as AODs and water vapor from direct- Sun measure- 
ments averaged for 16 minutes before and after the almu- 
cantar scan. As in Smirnov et al. [2003a], in this study the 
level 2,0 AERONET inversions are filtered to select only 
those retrievals which likely correspond to clean maritime 
aerosol. These constraints are based upon the measured 
AOD and retrieved size distribution. Here, t\ denotes the 
AOD at wavelength A (in nm), and a the Angstrom param- 
eter, which describes the spectral variability of r. Fi'ora the 
observations of Angstrom [1929] that ^ ,dA ''“ (where the 
turbidity coefficient I3 is the AOD at 1 pm), typically the 
definition a = -dln(r)/dln(A) is used, and a calculated 
between a pair of wavelengths Ai, A 2 as 


[ln(TAi) - 1“(^2)] 
[In(Ai) - ln(A 2 )l ’ 


(8) 


which, as a ratio of logarithms, is independent of logarith- 
mic base. The parameter Q' 44 (), 87 o is a standard AERONET 
product, and from this point all references to a indicate 
fi 440 ,a 70 - In the AERONET record a least-squares fit of all 
AODs within that spectral range is used to calculate a to 
reduce the impact of noise; otherwise, at low optical depths 
the uncertainty on direct-Sun r (of order 0.01; Holhen et al..^ 
1998, Eck et al., 1999) can propagate into significant un- 
certainties on a (see examples given by Wagner and Silva, 
2008). 

The constraints are that rsoo < 0.2 (where, if not avail- 
able directly, T 500 is estimated from the nearest available 
AERONET wavelength and a using Equation 8), and that 
0-1 Y o? < 1. These eliminate cases where there is likely 
residual cirrus contamination or some non-maritirne com- 
ponent in the aerosol, such as a local or transported pollu- 
tant, as pollutants are typically fine aerosol particles with 
large positive a, while desert dust and cirrus clouds have 
small or negative a, and the background maritime AOD 
is typically low (for example Eck et al, 1999, Ka.ufman 
et al, 2001, Dubovik et al, 2002, Knobelspzesse et al, 2004, 
Smirnov et al., 2009, 2011). Smirnov et al. [2003a] used 
^500 E 0*15 and a < 1; in this work, the rationale be- 
hind altering the constraints was bavSed on manual inspec- 
tion of size distributions which revealed that cases where 
0.15 < Tsoo < 0.2 generally still appeared maritime in 
character (and increased the data volume by approximately 
10%), while the lo'wer bound on a removed a small number 
of distributions which appeared suspect. However, if the 
thresholds used by Smirnov et al [2003a] are retained, re- 
sults are numerically very similar. The remaining retrieved 
size distributions are classified as ‘pure maritimef and the 
number of such distributions, f.ts well as the proportion of 
all distributions meeting these criteria, are given in Table 1. 

Midway Island has, of the well-sampled sites, the highest 
proportion (86%) of retrievals meeting the maritime crite- 
ria, and Kaashidhoo the lowest (20%). The average .size 
distributions for all retrievals, and all retrievals designated 
pure maritime, calculated in each case from the median 


d V^(r)/din(r) in each size bin, are shown in Figure 2 for 
these sites. The difference is minimal at Midwmy Island, 
which is not surprising. At Kaashidhoo, the coar.se mode 
is not significantly different between the two cases. How- 
ever, the ‘all data' average shows a significantly stronger 
fine mode contribution (with a peak around 0.2 jam), corre- 
sponding to transported continental aerosol. For both sites, 
the range of dV{r)/d\n{r) is large, particularly near the 
fine-mode and coarse-mode peaks. 

It should be emphasised that these size distributions still 
represent a retrieval, rather than direct measurement, of 
aerosol properties. Due to the selection of low-AOD cases, 
AERONET-retrieved refractive indices are not reliable in 
these situations {Dubovik et al, 2000) and so are not con- 
sidered. Nevertheless, AERONET offers a much greater sen- 
sitivity to aerosol parameters than current satellite instru- 
ments, and provides the most comprehensive ground-based 
dataset £tvailable, in terms of spatial and temporal coverage, 
data quality, and consistency of calibration and processing. 
Additionally, the large sample size, use of media.ns to de- 
crease sensitivity to outliers (from retrieval error or residual 
non-marine conditions), and fact that the inversions u.sed 
p^iss the AERONET quality control criteria to be raised to 
level 2 {Holhen et al, 2006), mean that the size distributions 
considered should be suitable for quantitative analysis. 

2»2. Properties of average size distributions 

For each of the 22 size bins, the median volume density 
from those inversions meeting the pure maritime criteria has 
been extracted to define an average size distribution for each 
site. This is what is meant by the term ‘median' or ‘aver- 
age' distribution through this work (i.e. median calculated 
for each individual size bin, rather than the median total 
aerosol volume). The use of medians is to limit the sensi- 
tivity to outliers. If means are used instead then the results 
are insignificantly affected at most sites (although the total 
aerosol volume typically increases, as most of the outliers 
are of higher- than- average volume). Throughout this w^ork, 
if a median is taken of a set containing an even number of 
values then the numerically larger is taken; this choice hjis 
a negligible impact on the results. The averaged distribu- 
tions are shown in Figure 3. The sites all show a similar bi- 
modal distribution, with a fine mode peaking at 0. 1-0.2 pm 
and a coarse mode peaking near 3 pm. Visually it resem- 
bles a bimodal lognormal distribution, although the coarse 
mode is persistently skewed, with a wider tail on the low- 
radius end. The broad similarity between sites is an indica- 
tion of the similar origins of the aerosol in different global 
oceans, although the abundance of both modes can vary, 
with Ascension Island and Midway Island having notably 
higher coarse-mode volumes than the other sites. Crozet 
Island has the smallest aerosol volume, perhaps because of 
the site’s elevation; it is also amongst the most skewed of 
distributions. 

The AERONET level 2.0 inversion product includes esti- 
mate es of Cv, i\, and (Xv (as well as rcff) for the overall size 
distribution, as well as, separately, fine and coarse modes 
(hereafter denoted by subscripts f and c respectively). How- 
ever, these calculations follow from Equations 3 and 4, with 
the .separation between fine and coarse modes determined by 
estimating the inflection point in the retrieved binned size 
distribution, fis opposed to from a fit to an assumed dis- 
tribution shape. The average of these parameters Is given 
in Table 2 for each site. The uncertainties are given as the 
scaled median absolute deviation (SMAD) about the median 
for each parameter, 


SMAD(x) = /3(x-x), (9) 

where ~ indicates a median quantity and A is a scaling fac- 
tor. If the underlying distribution is Gaussian, then the 
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SMAD is equivalent to the standard deviation for Q — 
1.4826, which is assumed here. The use of SMAD rather 
than standard deviation here is again to reduce the sensitiv- 
ity to outliers, and so provide a more representative estimate 
of the variability of the aerosol size distribution parameters 
in clean maritime conditions. 

Given the skewedness exhibited in Btgiire 3, the aver- 
ages of the size distribution parameters provided in the 
AERONET product are not the same as best-fit lognormal 
distribution parameters from the averaged size distribution. 
As a distribution constructed from lognormal components 
is desired, several approa.ches to addressing this tlierefore 
suggest themselves: 

1. Use the AERO NET- derived parameters from Table 2 
directly as btmodal lognormal distribution parameters, even 
though the underlying distribution is skewed. This will be 
referred to as the 'AERONET lognormal' method. 

2. Fit the average size distribution to a bimodal lognor- 
mal distribution. For this tnethod, bins up to and including 
the inversion point are lused to fit the fine mode, and bins 
with larger r the coarse mode, with a least-squares equal- 
weighting method. This will be referred to as the 'fitted 
lognormaf method. Although the bins adjacent to the inver- 
sion point may contain volume from both modes, in practice 
their inclusion or omission was found to have an insignificant 
effect on the fits. 

3. Fit a trimodal lognormal distribution to the retrieved 
AERONET size distribution, wTere the larger two modes 
will represent the skewed coarse mode of a bimodal distri- 
bution. 

The objective is then to deteirmine which of these methods 
leads to a distribution whose radiative properties match best 
those of the (non-lognormal) AERONET observations. The 
preferred method is the simplest one to match within the 
uncertainties of the data. Trimodal distributions are not in- 
vestigated in detail here; succeeding sections of this work mil 
show that a bimodal distribution is sufficient and the added 
complexity of a trimodal distribution is not required for this 
particular application (although the addition of modes does, 
of course, improve the fit to the retrieved size distributions). 
Ill the studies cited throughout this work, both bimodal and 
trimodal aerosol models have been used, although bimodal 
are the more common. An example of the average retrieved 
distribution, and multimodal lognormal distributions fit to 
it, is shown for Lanai in Figure 4. The SMAD of each bin is 
generally proportional to dV'(?')/dln(r), \vhich occurs due 
to changes in Cva and Cv,c fo^* diEerent AERONET obser- 
vations. This suggests the primary variability between dif- 
ferent size distributions is ‘up-down’ (total volume) rather 
than ‘left-right’ (peak position). LIsing the AERONET dis- 
tribution parameters directly for a lognormal distribution 
results in a larger-r fine mode and srnaller-r coarse mode 
peak than the averaged retrieved AERONET distribution. 
The fine and coarse mode volumes obtained from both meth- 
ods are very similar. A trimodal lognormal fit is also sho%vn, 
which provides a very close match to the average distribu- 
tion. 

The lower portion of Table 2 show's the volume size dis- 
tribution parameters for each AERONET site for a bimodal 
lognormal fit, along v/ith uncertainty estimates. For both 
these and the measured parameters provided by AERONET 
(Table 2). the parameters are close for most sites, particu- 
larly those w/ith a high number of AERONET inversions, 
suggesting that an approach to define a global maritime 
fme-mode and coarse-mode may be successful. DiEerences 
between sites may arise from factors such as diEerences in 
typical wind speeds or humidity between the sites (discussed 
below). 

The largest coarse mode radii are found at Nauru. Hen- 
derson et al. [2006] report that wind-induced v/ave- breaking 
leads to sea salt aerosol production in the coastal surf zone, 


and formation of downwind cloud trails, at this site. It is 
possible that the larger radii observed by AERONET here 
arise as a characteristic of this wave- breaking, or the fre- 
quent cloud trails mean residual cloud contamination is more 
likely. Henderson et al. [2006] also note that wind at Nauru 
is predominantly from the east; it is possible that an asym- 
metric aerosol field could lead to a bias in the AERONET 
inversion. Because of this strong surf zone source, results at 
Nauru may be less representative of the open ocean. 

The fine and coarse eEective radii are 0.132 /mi and 1.70 
/im for the ‘AERONET lognormal’ method (i.e. Equation 
7 applied directly to AERONET distribution parameters), 
ancl 0.122 (im and 1.80 /im for the bimodal lognormal fit, 
as compared to 0.150 /am and 1.87 /im when calculated 
from the AERONET distributions directly (i.e. Equation 5). 
The ‘AERONET lognormaT method is therefore closer for 
the fine-mode eEective radius, but the bimodal lognormal 
fit is closer for the coarse mode. Overall, these distribu- 
tion parameters (for both methods) are within the range 
of other studies (such as summarised by Silva et al, 2002, 
Smiinov et al, 2002, or Smirnov et al, 2003a). As also 
noted by Smirnov et at [2003a] and Ahmad et ai [2010], the 
AERONET size distributions are narrower than the older 
models presented by Shettle a,nd Fenn [1979], 

Table 2 show's both mean and number- weighted- mean 
size distribution parameters. However, from this point the 
number- weighted values will be used, such that the influence 
of poorly-sampled sites which more frequently report outly- 
ing values (Kaashidhoo, Amsterdam Island, and Crozet Is- 
land) is mitigated. The same conclusions broadly hold if the 
unweighted multi-site mean is used instead, as the weighted 
and unweighted means are similar. 


The results in Table 2 are averages for all conditions. 
However, it is known (discussed later) that the aerosol load- 
ing is infiuenced by meteorological factors such as the wind 
speed or availability of moisture. In this section, the data 
are therefore examined for these relationships. The Na- 
tional Center for Environmental Prediction (NCEP) reanal- 
yseSj available at U horizontal resolution output every 6 
hours, are used for this purpose (Derber et al, 1991). Al- 
though coarse-resoliition, such mode! output products are 
nevertheless the only current source of global spatially and 
temporally complete meteorological data, and so the only 
recourse if such information is to be used as an input to a 
global multiyear satellite aerosol retrieval scheme. 

Wallcraff, et ai [2009] compared satellite, buoy, and nu- 
merical w^eather prediction (NWP) model (including NCEP) 
near-surface wind speeds and found each dataset exhibited 
the same spatial patterns, although with regional relative 
biases. NCEP was found to have a positive bias relative 
to QuikScat satellite data, of up to 2.5 ms“^ in high-wind 
regions in the storm tracks. However, the correlation with 
QuikScat was high (greater than 0.9 over most of the open 
ocean) and bias smaller than 0.5 ms“^ in many oceanic re- 
gions, including those where most of the sites used in this 
study are located. The bias and root-mean square diEer- 
ence betv/een NCEP data and buoys were found to be small 
(0,15 ms~^ and 0.97 ms~^ respectively). The overall quality 
of agreement was similar for the diEerent datasets assessed. 
It is therefore suggested that any of these NWP datasets 
would provide a reasonable assessment of the near-surface 
oceanic wind speed near these sites. Wallcraft et al. [2009] 
also note that NWP results are likely to improve for more 
recent years than the period analysed (1998-2002) due to 
increased availability of observational data for c^similation. 


3. Relationship with meteorological parameters 

3.1. Data source 
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Figure 5 shows the results of comparisons between lin- 
early spatio-temporally interpolated NCEP wind speed and 
relative humidity near-surface fields with meteorological 
data recorded at approximately 10-minute intervals from 
twelve cruises of the Research Vessel (RV) Polarstern. These 
cruises are Atlantic Ocean transects, chosen to cover the lat- 
itude range inhabited by the relevant AERONET sites, and 
sampling a wide range of wind speeds [Konig-Langlo, 2005a, 
2005b, 2005c, 2005d, 2006, 2008a, 2008b, 2008c, 2009a, 
2009b, 2010a, 2010b). Overall the correspondence is high, 
although the coarser spatial and temporal resolution of the 
NCEP data mean there is a tendency for extrema to be 
missed. When the instantaneous wind speeds are averaged 
to daily values, the gradient of the least-squares best fit line 
forced through zero does not change much (0.84 to 0.S5), 
i.e. the NCEP data tend to underestimate the wind speed. 
However, the correlation increases from 0.90 to 0.96. For 
relative humidity, the data are almost unbiased (gradient 

0. 97) although the correlation coefficient is lower (0.54), re- 
flecting the typically small variability of relative humidity 
observed. These results support the validity of the use of 
NCEP data for the analysis of the relationship of aerosol 
properties with meteorology. However, the differences un- 
derscore the fact that analyses of this type are sensitive to 
not only the quality of the aerosol data, but also the mete- 
orological data. Part of the discrepancy may be due to the 
altitude difference (10 in for NCEP, as compared to 25 m 
above sea level for the ship). 

Additional AERONET aerosol products may provide fur- 
tlier insight into the relationship between marine aerosol and 
the ambient conditions. Firstly, although the AERONET 
size distribution inversions include temporal averages of 
direct-Sun AODs (as discussed previously), the full direct- 
Sun dataset is significantly larger. This is because direct- 
Sun measurements are taken every 15 minutes, while the al- 
mucantar scan is performed hourly, and requires completely 
cloud-free skies in the scanning plane. The second avenue 
is through the spectral deconvolution algorithm (SDA) data 
product, which provides the partition of AOD at 500 nm 
into separate contributions from the fine and coarse mode, 
and is independent of the other AERONET aerosol retrieval 
algorithms ( O^Neill et aL, 2003, 2006 for the current version 
4 level 2.0 dataset). 

Both of these additional products are therefore consid- 
ered. As before, the restrictions that rsoo ^ 0.2, 0.1 < a < 

1, and data from 1999 onwards (2010 for Graciosa), are im- 
posed. Additionally, to decrease the noise, and because of 
the coarser NCEP resolution, after obtaining the meteoro- 
logical information for each case, the direct-Sun and SDA 
products are downsampled to daily avei'ages for the com- 
parisons with wind speed and relative humidity before these 
thresholds are applied. This provides between 125 (Kaashid- 
hoo) and 1,171 (Midway Island) dates per site, with typi- 
cally around 15 measurements contributing to each daily 
average. This daily averaging is not done for the size dis- 
tribution inversions, as these are less frequent (many days 
have multiple direct-Sun meirisurements but no size distri- 
bution inversions). Additionally, the data are not averaged 
for the Y/ater vapor comparisons, as these are part of the 
AERONET product and so do not require matching with 
additional NCEP data. 

3.2. Dependence on wind speed 

Studies from both satellite and ground-based data have 
shov/n that increased near-surface Vvind speeds (u;s) are as- 
sociated with an increase in AOD, due to wind increasing 
the flux of sea spray (i.e. increased aerosol mass), and wa- 
ter vapor (i.e, aerosol swelling), from the ocean to the ma- 
rine boundary layer (recently Zidmski and Zielmski, 2002, 


Smirnov et aL, 2003b, Satheesh et aL, 2006, Mulcahy et aL, 
2008, Sakerin et aL, 2008, Glantz et aL, 2009, Huang et aL, 
2010, Lehahn et al, 2010, Adames et aL, 2011, Gramdey 
et aL, 2011, Kiliyanpilakkii and Meskhidze, 2011). These are 
similar to earlier results which directly sampled aerosol par- 
ticles, rather than remotely-sensed AOD (e.g. Lovett, 1978, 
Blanchard and Woodcock, 1980, Monahan et ah, 1983, Ex- 
ton et aL, 1985, Hopped et aL, 1990). Stronger wind speeds 
will, however, mean the aerosol is transported downwind of 
its source more rapidly, and so simple relationships between 
the two quantities are imlikedy to capture most of the vari- 
ability in aerosol loading, unless meteorological conditions 
are homogeneous over a large area and time period, and 
there are no other aerosol sources. 

Stronger correlations have been observed using the wind 
speed averaged over some time before the aerosol measure- 
ments were made (‘wind speed history'), rather than in- 
stantaneous wind sx:>eed, as wind may change on timescales 
shorter than aerosol lifetimes: the strongest correlations are 
typically found with wind speed averaged for 12-24 hours 
prior to the AOD measurement {Gathrnan, 1983, Hoppel 
et aL, 1990, Smirnov et aL, 2003b, Lehahn et aL, 2010). 
Some of these analyses bin data by wind speed, and then fit 
binned averages; this binning will naturally lead to stronger 
correlations than fits using all data points, as the variability 
is somewhat averaged-out, so this should be borne in mind 
when exarniiiirig regression statistics from dilferent studies. 
Mulcahy et al. [2008] found, for stable wind conditions, a 
very strong relationship between bin-average AOD and the 
approximate square of the wind speed. However, this was 
based upon measurements at a coastal site, and it is pos- 
sible that breaking waves on the rocky shore would lead to 
the production of additional aerosol above what w'ouid be 
observed in open-ocean (i.e. rock-free) conditions, or that 
there are differences in wind conditions beUveen the coast 
and nearby ocean (Blanchard and Wdodcock, 1980, Hen- 
derson et aL, 2006). It is uncertain whether these trends 
continue for high wind speeds, due to a paucity of data for 
W8 « 10 ms”^ or higher, and the few observations in these 
conditions have shown either increases, levelling-off, or de- 
creases in aerosol loading in strong winds (Blanchard and 
Woodcock, 1980, Exton et aL, 1985, Mulcahy et al., 2008, 
Pant et ad, 2008, Grandey et aL, 2011, Kiliyanpilakkii and 
Meskhidze, 2011). 

In this analysis, for each AERONET size distribution, 
NCEP near-surface (10 meter) wind speed data points are 
linearly interpolated in space and time to provide the 4n- 
stantaneous’ wind speed. The wind speed history is then 
defined by I'epeating this procedure backwards in time in 6 
hour increments, and averaging the resulting wind speed for 
up to 96 hours prior to the time of the AERONET retrieval. 
Pearson’s linear correlation coefficient between aerosol pa- 
rameters and this time-averaged wind speed are shown in 
Figure 6, separately calculated using all AERONET data 
together, and then individually for those sites with 100 or 
more AERONET inversions. The correlations v/ith the size 
distribution x>arameters are calculating using those directly 
reported by AERONET (i.e. those from which the avera,ges 
in the upper part of Table 2 were computed) rather than 
those from lognormal fits; similar results are obtained if the 
lognormal fit results are used instead. 

The strengths of correlations, and the time period which 
provides the strongest correlation, varies between sites. 
However, generally all sites show the same sign of correla- 
tion between the aerosol parameter and wind .speed history. 
Even if the true response of aerosol to changes in wind speed 
were the same at each AERONET site, the strengths of cor- 
relations might be expected to vary due to factors such as 
how well the coarse-resolution NCEP data represent the real 
wind speed: the range of wind speeds observed at a given 
location (if the response of the aerosol is nonlinear); and 
the fact that the sites are above sea level, so if a significant 
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proportion of the aerosol response takes place near the air- 
sea interface below the site then this may be missed in the 
AERO NET data if the aerosol is not vertically well-mixed. 
Blanchard and Woodcock [1980] present a model for the ver- 
tical dependence of sea salt concentration, based on wind 
speed; according to this, the highest number of particles are 
found at heights of up to 0.2 m above the sea surface, but 
most of these particles rapidly fall back in. LMr lieights of 
meters to several hundreds of meters, there is a small decline 
with height, with an eventual inversion layer of increasing 
concentration in the range 300 m to 600 m. Non-sea-.salt 
components were not considered. It is therefore likely that 
altitude and background aerosol contribute to the differences 
between sites. The issue is complicated by the fact that the 
results of Blanchard and Woodcock [1980] were collected in 
cloudy conditions, while the AERO NET data are clear-sky: 
Blanchard and Woodcock [1980] note that the salt inversion 
layer may depend on boundary- layer cloudiness. 

Positive correlations are found between the wind speed 
history and aerosol volume, particularly for the coarse mode, 
consistent with previously-mentioned studies. This mani- 
fests in additional positive correlations with AOD, stronger 
for tio 2 o fhan T 44 o, and a negative correlation with a, all 
linked to the fact that the fine mode is more optic ally- active 
in the visible, and the coarse mode in the near-infrared. 
There are also positive correlations between fv and the wind 
speed history, while the correlation is positive for the fine 
mode spread but negative for the coarse mode. These corre- 
lations are, however, generally weaker than those observed 
for Cv, r, and a, which themselves are typically 0.4 or 
smaller. 

Given the observed correlations, the next step is to exam- 
ine the size of the response of the aerosol size distribution to 
wind speed. For this purpose, averaged size distributions (as 
described previously) have been calculated by binning the 
AERO NET inversions according to the NCEP wind speed 
rather than by site. Bins have been chosen such that a large 
number of inversions fall within each, although there were 
only 67 cases of winds stronger than 10 ms '^S meaning care 
should be taken when considering results for high winds. 
Removing the constraints on rsoo and a at Lanai and Mid- 
way (the most well-sampled sites) does not result in signifi- 
cantly more high-wind points, implying that these imposed 
constraints are not causing a bias to low wind speeds. Ap- 
proximately 83% of inversions were for wind speeds from 
4-10 ms''h The resulting size distributions are shown in 
Figure 7. Table 3 details the number in each bin, r, and o, 
and shows the expected increase of AOD with wind speed 
and corresponding decrease of o. Size distribution parame- 
ters (for both the median of the corresponding AERONET 
inversion parameters, and lognormal fits to the median of 
distributions) are given in Table 4. The highest winds have 
a slightly lower AOD than anticipated by this trend, due to 
a slight decrease in the fine mode volume, although as men- 
tioned sampling is comparatively poor for this range. The 
base AOD for the calmest waters appears to be 0,068 at 440 
nm and 0,04 at 1020 nm. This is similar to (but slightly 
higher than) background AOD at 870 nm for dust- free pe- 
riod at Barbados of 0.035-0.04 reported by Smirnov et ai 
[2000b], The results for typical wind speeds also match well 
the ‘baseline maritime’ AOD at 500 nm of order 0.052-0.071 
reported by Kaufman et ai [2001], if the values in Table 3 
are interpolated spectrally, and observations taken on the 
decks of ships {Smirnov et al.. 2011). 

It is observed that, as wind speed increases, fine mode 
properties show mixed trends. The change in Cv,f is small, 
but potentially a decrease as wind speeds increase. The 
coarse mode exhibits a larger change: Cv,c increases strongly 
with higher winds, consistent with the previously-mentioned 
studies- The difference between bin averages for the lov/est 
and highest winds is a factor of two. Additionally, ry,c in- 
creases, although (7c varies less strongly. The increase of 
reff,c is driven mostly by the changes in ry,c. 


The difference between Tv.c for the two most populated 
bins (4-6 rns“" and 6-8 ms~^) is within the variability of 
distributions in each bin (for the ‘AERONET logriormaL 
method) and smaller than the fit uncertainty on each bin 
(for the dognorraal fitted’ method). The same is true for dc. 
They are also smaller than or comparable to the variability 
or fit uncertainty of these parameters at individual sites (Ta- 
ble 2). Additionally, the values of these parameters for the 
0-4 ms~^ and 8-10 ms“’^ bins are typically within or close 
to the variability or fit uncertainty. These are important re- 
sults as they indicate that file multi-site average -rv.c and ctc 
may be able to represent coarse-mode aerosol for the major- 
ity of encountered wind strengths. This highlights again the 
underlying similarity of the coarse-mode aerosol at differ- 
ent locations (i.e. aerosol found at different locations with 
.similar wind spetxls correspond more closely than aerosol 
observed at a single location over a range of wind speeds). 
The change in the results of this analysis are small if size 
distributions are binned by the 24- hour- average, rather than 
instantaneous, wind speed. 

3.2.1. Fits of Cv to wind speed 

Linear and exponential fits of AERONET retrieved fine 
and coarse mode volumes (Cy,f and Cy.c) to the 24-hour- 
averaged wind speed are shown in Figure 8, for data from 
all sites considered together. In both crises, the correlation is 
much stronger for the coarse mode than fine, although both 
linear and exponential models provide a similar quality of 
fit, due to a comparative lack of data for high wind speeds 
where nonlinearity would be more evident in the exponen- 
tial model For the most commonly-encountered wind speed 
ranges, both methods give very similar results. The fine 
mode volume is iridependerff of wind speed, while the coarse 
shows a positive correlation, consistent with the mechanism 
of wind-driven emission. The 24-hour-average wind speed 
is used as it shows a stronger correlation with aerosol vol- 
ume than instantaneous wind speed in most cases (Figure 
6), although coefficients of fit are similar if instantaneous 
wind speed is used instead. Similar results are obtained if 
the fits me performed on a site-by-site basis (omitted for 
brevity), with the fine mode independent of wind speed and 
the coarse mode typically wuth a base volume of 0.01-0.02 
and an increase of order 0.005 ^m'^/nn '" per 1 
ms'“^ wind speed, although the smaller individual sample 
sizes lead to high uncertainty on fit coefficients. 

For the fine mode in Figure 8, high outliers contribute 
to the fact that the least-squares regression linear inter- 
cept (0.0064 /im'“^) is higher than the average observed 

fine mode volume for a lognormal distribution (0.0057 

or 0.0056 /:xm^ from upper and lower portions 

of Table 2). For these very low Cvj, least squares regres- 
sion is less appropriate as the spread of values may not 
be Gaussian (i.e. the AERONET data will not retrieve a 
negative volume; thus, noise on the low-volume retrievals 
will be on average biased positive). If an alternative fitting 
method than least-squares is used, such as least absolute 
deviation (Bloomfield and Steiger ^ 1980) or reduced major 
axis {Ayers, 2001), then the fine-mode fit parameters change 
(although correlations remain effectively zero) while coarse 
mode fit parameters are affected minimally in most cases. 

Averaging the data over all sites in bins of 0.5 ms“^ pro- 
vides significantly stronger correlations, as shown in Fig- 
ure 9, for linear relationships betw^een wind speed and Cv 
(equations of fit shown in the Figure). To reduce the effect 
of outliers, bin medians rather than means are used. Sev- 
eral bins at very low and high winds contained very small 
numbers of retrievals, meaning the bin averages were noisy, 
a.nd as these bins lie at the edges of the range they have a 
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strong influence on the linear least-squares fit. For this rea- 
son, bins with fewer than ten size distributions have been 
excluded from the analysis in Figure 9. For the fine mode, 
ciuantisation in the AERONET Cv (increments of 0.001) is 
evident in the fit. Bt^tween the quantisation and variability 
within each bin, there is effective no relationship for the fine 
mode. 

These results imply that, if only the wind speed or wind 
speed history is known, fussuming ^ 0.0056 (i.e. the 
average observed value) and Cy,c 5=^ 0.02 + 0.0025a’s or 
Cv,c ^ 0.015 + O.OQZCnjJS (i.e. the global average best lin- 
ear relationship for the binned data, dependent on w^hether 
last ant aneous or 24-hour-averaged wind speed is known) will 
give a reasonable first-order estimate of the aei'osoi volume. 
This is examined further in Section 5. Despite this, wind 
speed alone is likely to be a poor predictor of aerosol volume; 
the variability within bins on Figure 9 is similar to the range 
of volumes encompassed by the best- fit line. Further, the 
different gradients of coarse-mode volume with wind speed 
presented in this section vary by a factor of two, as a result 
of simply changes in data aggregation. The data do not per- 
mit analysis for very high wind speeds. To an extent this 
will be related to inadequacies in the NCEP data represent- 
ing the true wind speed history at each site, and errors in 
the AERONET size distribution retrievals. However, this 
highlights the necessity for complexity and consideration of 
the aerosol life cycle from emission to removal in modelling 
of the aerosol burden, as is performed by the current gener- 
ation of global models. 

3.2.2. Fits of AOD to wind speed 

Statistics of linear fits of direct- Sun AOD and a to wind 
speed (in both cases, from points averaged for each day) 
are presented in Table 5. There is considerable variety be- 
tween the sites, both in terms of strength of correlation and 
the fit parameters, which may in part reflect different local 
sources. In general, the strongest agTeement is found be- 
tween ws and a; stronger correlations are found with rio 2 o 
than T 44 Q^ These results can again be explained in terms of 
an increased coarse-mode presence at higher wind speeds. 
Due to the high scatter, the uncertainties on these linear fit 
parameters (not shown) are large. There is effectively no 
significant correlation between tsqo.c and wind speed. Al- 
though this is a surprising result, because the coarse mode 
optical depth at this wavelength is low (as compared to 1020 
nm, where almost all the AOD comes from the coarse mode), 
it is likely that any signal is masked by the uncertainty on 
T 500 .C and wind speed, or background variability. Another 
possibility is uncertainties in the fine/coarse partition in the 
SDA data. The relationship between Tsoo,/ and wind speed 
is similar to that of r 44 o; this can be explained by the fact 
that AOD at 440 nm is mostly determined by the fine mode, 
so T44G and T 5 oo,/ are sensitive to the same parts of the un- 
derlying aerosol burden. 

Averaging the data over all sites in bins of 0.5 ms 'E and 
taking the bin medians, leads to the relationships .shown 
in Figure 10. Again, poorly-sampled bins (fewer than 10 
points) are omitted. An approximate linear relationship ap- 
pears to hold for all cases, although the largest outliers on 
the fit occur for wind speeds higher than 12 rns“^ As in 
the case for aerosol volume, the variability within each bin is 
similar to or larger than the range of the bin-average values, 
again illustrating that wind speed alone is of limited utility 
in predicting the marine aerosol burden for an individual 
case. 

3.2.3. Joint analysis with sea surface temperature 

Marine {sea, spray) aerosol source functions in Earth sys- 
tem models are typically the product of an assumed size dis- 
tribution with a (wind-speed-dependerit) fractionai white- 
cap cover, although there is considerable scatter between dif- 
ferent parainetrisations {Lewis and Schwartz. 2004, O'Dowd 
and de Leeuw. 2007). Jaegle et ai [2011] improved the corre- 
.spondence betvyeen modelled and observed sea salt aerosol 


burdens by modifying the source function of Gong [2003] 
with an empirical sea surface temperature (SST)-dependent 
correction. This dependence is thought to arise partially due 
to the changing kinematic viscosity of the sea surface with 
temperature (determining bubble rising and breaking), and 
may also partially be because observed whitecap fraction is 
also linked to SST {Anguelova and Webster. 2006). An SST 
dependence has also been observed in laboratory studies of 
seawater and analogues {Sellegri et al. 2006, and references 
therein). 

Due to diurnal changes in SST, aerosol lifetimes on the 
order of days, and the fact that the AERONET sites are 
not located at the ocean surface, the observed aerosol in the 
atmosphere at any given time may not be representative of 
the aerosol flux from the ocean for the temperature at that 
given time, and there is no direct match possible between 
SST and the AERONET inversions. For these ret\sons, the 
version 2 Optimal Interpolation (01) SST dataset {Reynolds 
et at, 2006) is used for a joint analysis of effect of wind 
speed and SST on aerosol. This provides global gap-filled, 
bias-corrected, daily average (daytime and nighttime orbits) 
bulk SST. As the SST is bias-corrected against buoys, this 
bulk SST corresponds to a depth of order 0.5 m below the 
surface, and is typically within 0.5 K of the skin SST, al- 
though this depends on meteorological factors (e.g. Murray 
et at.. 2000). The dataset is derived from Advanced Very 
High Resolution Radiometer (AVHRR) and Advanced Mi- 
crowave Scanning Radiometer (AMSR) data. It is provided 
on a 0.25^^ grid but here is downsampled to 2.5^ resolution 
to provide a better representation over the larger source re- 
gion that the AERONET site may sample from on a given 
day. 

The mean SST is 24.5'^C (median 24.7“C), and the stan- 
dard deviation 2.9"^C, The coolest and warmest tempera- 
tures encountered are 4.2'^C and fil.O'^C respectively, al- 
though the number of cases with water cooler than 20” C 
is very sniaiL This is because the majority of the sites are 
in warm tropical waters, and so any conclusions drawn may 
be unrepresentative of cooler waters. 

Figure 11 is analagous to Figure 7, except the data are 
also subdivided by SST, The 24- hour- averaged wind speed 
li£is been used to stratify the data, although the results do 
not change significantly if the instantaneous wind speed is 
used instead. The SST bins have been chosen to be nar- 
row while still retaining sufficient sampling in as many cases 
as possible, although this is difficult for the highest wind 
speeds {ws > 10 ms““^). Despite the previously-documented 
links between SST and marine aerosol production, there ap- 
pears no strong and consistent link with the size distribution 
here (certainly compared with the effects of wind speed). 
For high winds, the coolest (SST< 21”C) and w^armest 
(SST> 27” C) have a higher coarse-mode volume than the 
intermediate SST ranges, although the sample size for these 
wind and SST bins is very limited, so these results should 
be interpreted with caution. 

Given the rough linear relationship between wind speed 
and coarse-mode volume observed in the previous section, 
and the results of Jaegle et al. [2011] that scaling the source 
function for marine aerosol production by a cubic poly- 
nomial in SST improved the marine aerosol burden in a 
global model, similar fits (omitted for brevity) of the form 
Cv = (a + 6SST + c-SST^ + dSST'*) x (e + f x ws), where 
a — / are coefficients to be determined, have been attempted 
here, for fine and coarse mode volume. However, using such 
fits Pearsoirs linear correlation coefficient showed only slight 
improvement as compared to the wind-s peed-only relation- 
ship for both fine and coarse modes, strengthening from 0,03 
to 0.12 and 0.38 to 0.42 respectively. These results suggest 
that consideration of SST provides minimal additional skill 
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in. predicting marine aerosol loading. This may reflect un- 
certainties in how well the datasets (AERONET, wind, and 
SST) are able to represent the true aerosol or meteorological 
conditions, or simply that although SST may affect aerosol 
production, aerosol lifetime is sufficiently long that it has 
little impact on the total aerosol burden in these regions. 
The SST range sampled may be a significant factor, given 
that the largest corrections to the aerosol source function 
made by Jaegle et al. [2011] are outside the typical range of 
AERONET data used here. 

3.3. Dependence on relative humidity 

Historically, a common approach to modeling aerosol size 
distribution and refractive index, such as in the frequently- 
used database of Shettle and Fenn [1979) and those which 
draw from it (e.g. Hess et al, 1998), has been to initially 
define properties for a ‘dry’ aerosol type of some asvsumed 
composition. These ‘dry’ properties are then modified ac- 
cording to variations in relative humidity (rh) using models 
of aerosol swelling (e.g. Hand, 1976, Kotchenmther et al, 
1999 for marine cases). The effect of this is, as the relative 
humidity increases, the size distribution shifts to larger par- 
ticles and the refractive index approaches that of water, with 
the change dependent on composition. As the size distribu- 
tions in this work are calculated from AERONET inversions 
they represent the aerosol size distributions as found in the 
wild’, and their variability will encompass the effects of the 
range of relative humidity and consequent aerosol sw-elling 
and drying. There are difficulties when trying to use a 
methocl based on knowing the relative humidity to define 
an aerosol model. Analogously to wind speed, the quality 
at which coarse-resolutiori model relative humidity is able 
to represent the actual relative humidity on a finer scale 
is likely highly variable. Aerosol swelling shows hysteresis, 
such that even if the relative humidity were known this may 
not be sufficient to model the extent of swelling unless the 
prior history of the air mass is also known {Kotchenruiher 
et al, 1999). Additionally, there is evidence that in the case 
of sea salt aerosol the currently-used mixing rules do not re- 
produce the observed refractive index (Irshad et al, 2009), 
Finally, it should not necessarily be expected that the ADD 
should correlate well with near-st.irface humidity, as AOD is 
a columnar quantity. 

NCEP rh data are interpolated here to AERONET inver- 
sions in the same way as the wind speed data. There is lit- 
tle direct correspondence between the relative humidity and 
wind speed (not shown), and the interval 60 < rh < 80% 
contains 88 % of all relative humidities encountered. Figure 
12 shows the averaged AERONET size distribution (calcu- 
lated as previously) for inversions aggregated by relative hu- 
midity; Table 6 shows the number of retrievals in each bin, 
as well as r, a, and the mean wind speed for the data in that 
bin (which is similar for each, varying by 1.7 ms” ^ or less 
between bins). Table 7 gives statistics of these distributions 
and fits to them, in the same w^ay as Table 4 for wind speed. 

Some studies have observed an ariticorrelation of r and rh 
for rh <75%, and a positive correlation for rh >75%, when 
measured simultaneously on ships {Smirnov and Shtfrtn, 
1989, Srrmmov et al, 1995). This was attributed to tur- 
bulent exchange in the marine boundary increasing r and 
decreasing rh, leading to a natural anticorrelation, but this 
effect being overvvTielrned by moisture uptake for high hu- 
midities (with a deliquescence point at rh ss75%, depen- 
dent on temperature and composition). Upon first exami- 
nation, the combination of AERONET and NCEP data pro- 
vide some evidence to support this relationship. The AOD 
at 440 nra and 1020 nm, as well iis decrease with in- 
creasing rh until around rh -=60%, For rh >90?^ they (and 
Cvj) jump dramatically. Other size distribution parameters 


change less dramatically (except for the most extreme, and 
most poorly-sampled, bins). This behaviour is consistent 
with increased turbulence leading to increased particle num- 
ber (but little change in particle type) in ‘dry’ conditions, 
but rapid aerosol swelling in ‘moist’ conditions. There is an 
increase of a with rh up to rh :=85 %, and then a drop. The 
drop in a for the lowest and highest humidities could be re- 
lated to the swelling of aerosol but could also be an artefact 
of cloud contamination or limited sampling. The variability 
of O' within most bins is large. 

However, the evidence is weak, because of the poor sam- 
pling for low and high humidities. In particular there are 
only four cases of rh >90%, and the average size distribu- 
tion is unusual compared to the other ranges. The fact that 
these are grid-box average humidities suggests the presence 
of clouds in some region of the grid box is likely, and so it 
may be that these distributions are unphysical and contam- 
inated by clouds. All four of these cases occur for Graciosa, 
from which the data have not yet been raised to level 2.0 (i.e. 
full quality control has not yet been applied). Similarly, for 
the range 80 < r/i < 85%, 27 of the 38 points are from Gra- 
ciosa. Of the 71 points with rh <60%, 40 occur at Midway 
Island, 19 at Lanai, and 6 at Bermuda. Midv/ay Island and 
Bermuda can on occasion be influenced by transported dust 
{Srnim:OV et al, 2000b, 2003b) and so Lanai may also. It is 
therefore possible that the low-humidity results here are in- 
fluenced by dust transported in dry air masses, rather than 
a change in the abundance of marine aerosol All of these 
pointvS at Midway occur from December to April, when dust 
transport is expected to be most likely. If points from Mid- 
way and Lanai are removed, then for rh <60% the mean 
T 440 =0.068, no2o=0.047, and a'~0.38, although sampling 
becomes very poor. In these cases the AODs and .size distri- 
butions for low humidities match closely those for other lui- 
raidity ranges (except the suspect Graciosa rh >90% data) 
in Table 6, and the trend in AOD with rh is effectively 
removed, although the trend in a remains. Removing the 
Bermuda data does not have a significant effect on the re- 
sults. It is therefore possible that this small number of the 
driest cases represent residual contamination by transported 
dust. The coarse-rnode peak radius is also shifted to larger 
volumes for these drier bins, which supports this (Figure 
12). if these seasons are removed for these sites for the 
previous main analysis of aerosol size distributions (Section 
2.2), however, the impact is negligible. 

Some of the \^arlability in ail these cases will arise from the 
hysteresis of aerosol deliquessence (he. the path by which the 
current relative humidity v/as reached is important), which 
may mask any change in aerosol properties expected around 
rh =75%. Inadequacies in the coarse-resolution NCEP data 
will also limit the extraction of useful information. Over the 
heavily- populated range 60 < rh < 80% size distribution 
parameters show' little change, suggesting that average val- 
ues are sufficient to describe the majority of cases encoun- 
tered, The sampling is inadequate to make confident state- 
ments about behaviour in low or high humidities. To further 
examine these relationships between relative humidity and 
AOD for dry and moist air, linear regressions have been 
performed using the direct-Sun data. Binning daily average 
data from all sites together (as in Figure 10 for wind speed) 
leads to Figure 13, the trends of which are broadly con- 
sistent with Table 6, although again the variability within 
each bin is typically larger than the range spanned by bin- 
median values. Sampling is insufficient to examine data for 
rh >90%. The same conclusion is reached if other wave- 
lengtlis are used; whether the data are subset according to 
wind speed or not: whether instantaneous or daily-averaged 
data are used; or whether sites are considered individually 
or jointly. 

3.4. Dependence on water vapor content 

In addition to the aerosol bands, the Sun- photometers 
used in AERONET have a channel around 940 nra v/hich en- 
ables the retrieval of water vapor with an uncertainty of or- 
der 5%-10% {Smirnov et al, 2004, Alexa7idrov et al, 2009). 
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This provides an alternative way to examine the effect of 
moisture on maritime aerosol The columnar water vapor 
is provided in units of g cm~^ (equivalent to cm, given a 
density of 1 g cm"‘^). Separating the AERONET size distri- 
butions according to columnar water vapor gives average dis- 
tributions shown in Figure 14. The lowest bin (0 < wv < 1 
cm) has the low^est volume (and AOD), but aside from that 
there is no significant dependence of distribution parameters 
on water vapor content (tables omitted for brevity). Table 
8 shows the variation in AOD with columnar water vapor, 
along with the average wind speed for each bin; there is no 
trend. The case of the 0 < wv < 1 cm bin is sampled poorly, 
and approximately half (7/15) of the size distributions are 
from Crozet Island, so in this case the low AOD and water 
vapor may both be due to conditions specific to this site 
rather than an effect of water vapor on AOD. It is therefore 
likely that results for this water vapor range do not represent 
open ocean conditions well. There is a small increase of a 
with water vapor; however, the variability on a within each 
bin (0.15-0,2) is of similar size to the range over all bins. 

The relationship has also been examined for individual 
vSites, and restricted to different siibsets of wind speeds, to 
investigate whether aggregating data from multiple sites or 
wind regimes together was masking the signal. However, 
this did not reveal any significant relationships. The lack 
of correlation could be explained as a combination of ef- 
fects resulting from the low ranges of AOD and water vapor 
encountered; that the vertical distributions of aerosol and 
water vapor have small overlap thus limited potential for 
influence; or the possibility that the timescales of aerosol 
growth and water vapor variability are different. This site- 
dependent relationship (or lack thereof) between AOD and 
moisture has been noted in previous studies {Exton et al, 
1985, Hoppel et al. 1990, Smirnov et ai, 1995, 2000c, Holhen 
et al.y 2001, Sakerin et al, 2008). 

The strength of the linear correlation coefficient R be- 
tween water vapor content and AOD for direct-Sun data is 
0.2 or less in most cases when calculated for any site or 
wavelength, for a selection of assumed relationships (lin- 
ear, quadratic, or exponential). This confirms the results 
from the smaller AERONET inversion dataset that the in- 
fluence of water vapor on AOD in these pure maritime con- 
ditions is small In case the restrictions rsoo < 0.2 and 
0.1 ^ oc < 1 were masking the true relationship (by re- 
moving those points where elevated water vapor was asso- 
ciated with an increased AOD), results without these two 
constraints have also been examined. However, the rela- 
tionships remain weak; an example is shown for Lanai and 
Midway Island in Figure 15, for exponential fits between 
wv and Tf,rm or rc,boO‘ At Midway Island there is evi- 
dence that enhanced water vapor corresponds with a de- 
crease in AOD. This could be related to periodic transport 
of Asian dust in dry air masses (Smirnov et ai., 2003b, Eck 
et at., 2005) rather than an effect of water vapor on marine 
aerosol. Also shown in Figure 15 are analagous results for 
the AERONET sites of Wallops (coastal; 37.942*" N, 75.475*" 
W) and COVE (a platform 25 km from the coast; 36.900*" 
N, 75.710*" W). Because of their coastal (rather than remote 
island) locations they are more susceptible to continental 
influence, which is reflected in the higher AODs than ob- 
served at sites like Lanai or Midway Island. For both the 
‘ail points’ and ‘maritime conditions’ datasets, there is a sig- 
nificant strong exponential relationship between AOD and 
water vapor content, with R = 0.63 — 0.77 for fits to all 
points, and R — 0.2 — 0.44 for only those classified as mar- 
itime conditions, which is still stronger than most of the 
relationships found for the sites considered for this study 
(Table 1). This provides further evidence that the lack of 
correlation observed for the maritime sites under maritime 
conditions is real. 

3.5. Other factors 

Various studies have showm that the organic content of 
submicroo marine aerosol is increased in waters with high 


biological activity, signified by the enhancement of organic 
carbon within sea water (such as O’ Dowd and de Lccmuk 
2007, Fuentes et ai, 2010), although the relationships are 
complicated. Fuentes et at [2010] found, for experiments 
with seawater proxies enriched with algal species, that the 
number of generated aerosol particles of modal dry radius 
of approximately 0.02 jum was increased by up to approxi- 
mately a factor of two as compared to a proxy without the 
algae. The effect on particles of dry radius of order 0.05 pm 
and larger, which compose the bulk of the volume of the fine 
mode, and so fine mode optical depth, was smaller. There- 
fore it is unlikely that this enrichment has a strong effect 
on the visible AOD, Additionally, investigation of the effect 
on the size distribution using AERONET aerosol and satel- 
lite organic activity proxies (e.g, chlorophylLa) data is dif- 
ficult, for reasons including the physical separation between 
the Sun-photometers and the water, the possibility for con- 
founding effects from errors in atmospheric correction in the 
ocean colour products, and the heterogeneity and difficulty 
of retrieval of ocean colour products in coastal waters. For 
these reasons, the relationship is not examined here. 

Some other factors influencing marine aerosol production 
are disc\,issed by Levm and Schwartz [2004], but are either 
difficult to assess using available data, or likely to have a 
minor influence on the AOD, and so are not further con- 
sidered here. These include atmospheric stability, precipi- 
taition, surface- active materials (such as the aforementioned 
organic carbon), wave state, boundary layer height, fetch, 
salinity, and bottom topography. 

4, Refractive index and calculated AOD 

As well as the size distribution, knowdedge of the com- 
plex refractive index m ~ n~~ik, where n is the real compo- 
nent and k the absorption coefficient, is required to calculate 
the AOD at given wavelength. Although provided in the 
AERONET record, for the low AODs considered here the 
uncertainty on this is large and so results may be unreliable 
{Dubovik et al. 2000). Additionally, the inversion provides 
a single refractive index for the aerosol model as a whole. 
Observational evidence suggests that the fine mode is com- 
posed predominantly of a mixture of sulphates, organic com- 
pounds, and salt, while the coarse mode is predominantly 
salt, although the exact composition is variable and depen- 
dent on meteorological and biological factors (Hegg et ai, 
1997, Magi et ai, 2005, Clarke et ai, 2006, O’Dowd and 
de Leeuw, 2007, Fuentes et al, 2010), These different com- 
positions would be expected to lead to difference refractive 
indices for the two modes. For these reasons, a variety of 
refractive indices are tested in this work, shown in Table 
9, which includes ground-based observations as well as one 
pair of components used in the current Moderate Resolution 
Imaging Spectroradiorneter (MODIS) aerosol retrievM over 
ocean [Remer et ai, 2006), and the Optical Properties of 
Aerosol and Clouds (OPAC) database (Hess et ai, 1998), 
which is drawn from Sheitle and Fenn [1979], and used in 
various other satellite and model datasets (e.g. Sayer et ai, 
2010). This is by no means an exhaustive list, although 
it does encompass the range of common ly-used values; for 
example, Fraser et al. [1997] use rn = 1.38 — O.OOli, in- 
termediate between Smirnov et ai [2003a] and Silva et ai 
[2002] in Table 9. Although the refractive index has spec- 
tral behaviour, over the range of intere.st here most studies 
in Table 9 use a broadband value, as variation through the 
visible and near-IR is small. 

The refractive indices and size distribution fitting method 
are tested by attempting to recreate each case of AERONET 
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spectral AOD using Mie theory and the multi-site aver- 
age size distribution parameters rv,f, rv,c; and Cc from 
both fitting methods (Table 2). The distribution volumes 
Cv.f and for each case are taken from the AERONET- 
reported parameters (or a lognormal fit to the size distri- 
bution, as described previously) for ea,ch individual obser- 
vation. In this way the tests mimics the v/ay the average 
model may be implemented in satellite retrieval schemes, 
i.e. the spectral AOD is determined only by altering the vol- 
umes of each component while the distribution peak radius 
and spread are held constant. This allows an assessment of 
the degree to which the average model is able to represent 
maritime aerosol at each site, and will inherently include 
the effects of changes in wind speed, relative humidity, and 
errors resulting from potential aerosol nonsphericity. 

For each site, the correlation coefficient, median bias (cal- 
culated - AERONET observed AOD), and SMAD are cal- 
culated. The evaluation of each combination of size distri- 
bution and refractive index, how^ever, is restricted to only 
five sites in Table 1 with at least 100 observations meet- 
ing the maritime criteria (Lanai, Bermuda, Midway Island, 
Ascension Island, and Tahiti), as these each provide a repre- 
sentative dataset of maritime aerosol data of reasonable size. 
Although Gradosa and Nauru also have over 100 observa- 
tions, the former is omitted due to potential concerns about 
cloud-contamination, and the latter additionally due to sus- 
pected influence of surf-generated aerosol {Henderson et at. 
2006). Over this subset of sites, the minimum, maximum, 
and median of each of these parameters is presented for each 
case in Figure 16, This provides a simple reference of how 
well each potential combination of size distribution and re- 
fractive index is able to represent the AERONET AOD. 

Figure 16 reveals that, in general, the spread of statistics 
between sites is larger than the spread induced by changes in 
refractive index. All models tends to overestimate a, due to 
r 440 being comparatively unbiased while the AOD at longer 
wavelengtfis has a slight low' Him. For the purpose of satel- 
lite AOD retrievals, this is not likely to be a problem as the 
bi^is could simply be redressed by altering the total volume 
of the fine and/or coarse modes. However, this would mean 
that the bias w^ould translate from AOD into volume, which 
would then mean that derived aerosol mass estimates may 
be inaccurate. 

The ‘fitted lognormaL approach results in higher corre- 
lations between calculated and AERONET' AOD, with a 
lower spread of difference (SMAD). The correlation coeffi- 
cients are high in all cases, particularly considering that the 
range of AOD is smtill (most data are for 0.03 < r < 0.08). 
In contrast, this method leads to slightly more negative bi- 
ases in AOD. These biases typically remain smaller than 
the SMAD, and importantly both of these are often around 
the reported uncertaintv on AERONET AOD of 0.01-0.02 
(Holben et al, 1998, Elk et al, 1999). The SMAD is the 
metric of most interest because it provides information on 
the scatter of the AOD about this bias. For these reasons, 
the ‘fitted lognormaL method is deemed the more useful of 
the two approaches. This is an important result because 
it demonstrates that, just by varying the volume of each 
component, a single value of each of rv,!, rv,a <Tf, and <Jc 
are able to rep)roduce the AERONET AOD over a variety of 
oceans, wind speeds, and humidities, wdth a precision similar 
to that of the AERONET AOD measurements themselves, 
and even considering the fact there may be a non-maritime 
contribution to the aerosol loading in some cases. 

There is no clear ‘best case’ of refractive indices to choose. 
An examination of Figure 16 reveals that all cases tend to 
underestimate AOD from 675 nm onwards; for this wave- 
length, the majority of the calculated AOD conies from the 
coarse mode, while r is fine-dominated at 440 nm, and both 
modes contribute approximately equally at 500 nm (calcu- 
lations not shown). The similarity of the AOD bi<?is for all 
refractive indices for the longer wavelengths reveals that 


coarse-mode AOD is (over the range of typical refractive 
indices listed in Table 9) comparatively invariant with re- 
fractive index, which implies that choice of refractive index 
is probably not significant for background coarse- mode ma- 
rine aerosol. Conversely, for the fine mode those cases with 
a larger refractive index (2, 4, and 5, using Remer et al. 
2006 and Shettle and Fenn, 1979/ Hess et al., 1998) result in 
a higher AOD (and so less negative bias) at 440 nm and 500 
nm. However, this means they provide a larger overestimate 
of a. Given the low AODs eticountered, this overestimate 
of a is not considered problematic, as large errors in a can 
propagate from smaller errors in AOD {Wagner and Stlva, 
2008). Unfortunately, the AOD is not measured at longer 
wavelengths at these sites, wdiich means fhe applicability of 
the model (and any spectral behaviour of refractive index) 
at other common satellite wavelengths, such as 1.6 f.itn and 
2.1 pm, may not be assessed directly. 

Based on these factors, from this point case 4 from Ta- 
ble 9 (fine mode rn = 1.415 - 0.002L coanse mode ?n = 
1.363 3 X 10”^'^i) is used, although results change mini- 

mally if case 5 is used instead, or case 2 (for wavelengths 
of 675 nm or longer). The correlation, median (calculated 
- AERONET) bias, and SMAD for each site for this choice 
of refractive index, and the recommended ‘lognormal fitted’ 
distribution approach, are provided for r\ and a in Table 
10. These refractive indices (with size distribution compo- 
nent parameters fvj = 0.157 pm, ?v.c ~ 2.58 pm, (jf — 0,50, 
(7f — 0.72) are hereafter referred to as the ‘recommended 
maritime model’. The single scatter albedo is approximately 
0.98 over this wavelength range. 

The site with the most negative bias in AOD (for all re- 
fractive index cases) is Kaashidhoo; however, as discussed 
previously, this site is particularly poorly-sampled, and sub- 
ject to potential seiisonal influences of aerosol outflow^ from 
the Indian subcontinent, and so likely less representative 
of clean maritime conditions. The next-largest biases are 
for Graciosa, which may also potentially be influenced by 
transported continental or local pollution, or cloud contam- 
ination, as discussed previously. For both of these sites it 
is likely that a pollutant would be more absorbing than the 
background maritime aerosol, such that the maritime model 
would underestimate the fine-mode absorption AOD, wliich 
is consistent with the observed underestimates. Crozet Is- 
land has the largest scatter, although with a sample size of 
only eight, is very susceptible to outliers. The low AODs and 
high scatter at Crozet Island mean that a is reconstructed 
poorly. 

The calculations have also been performed for AOD at 
340 nrn and 380 nm (and are included in Table 10) which, 
although not used for the AERONET inversion, are available 
for parts of the record. The results for these wavelengths are 
omitted from Figure 16 as they are very similar to those 
at 440 nm, with spectral trends continued (e.g. stronger 
negative biases at Kaashidhoo). No additional insight into 
the most suitable refractive index is obtained. Additionally, 
when shorter wavelengths are considered the assumption of 
a spectrall^^-neiitral refractive index is also likely less appro- 
priate. 

5. Prediction of Maritime Aerosol Network 

AODs 

The previous sections have focussed on results obtained 
for coastal island sit€?s. The Maritime Aerosol Network 
(Mx\N) component of AERONET {Smirnov et al.. 2009) 
provides spectral AOD measured on ocean cruises using 
hand-held Microtops II Sun-photometers, with an uncer- 
tainty of approximately 0.02 (Knobelspiesse et aL. 2004). 
The AOD is reported at a selection of wavelengths (440, 
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500, 675, and 870 rim for the cruises used here), Using 
these data helps to establish the utility of the wind-speed- 
dependent relationships observed in previous sections. For 
this analysis, level 2,0 (cloud-screened and quality-assured) 
AODs from the ‘series- average’ product are used. One mea- 
surement series is defined in this product as the set of AOD 
measurements taken with a gap of no more than 2 minutes 
between an individual pair. 

A subset of cruises whose measurements took place in ar- 
eas likely to have minimal influence from transported aerosol 
sources are used. These are the SA Agulhas (during 2007- 
2008), Knorr (during 2008), Marion Dufrene (a cruise from 
each of 2008, 2009, and 2010), Sonne (during 2009), and 
Astrolabe (during 2009-2010). The MAN data are subject 
to the Tooo < 0.2 criterion to improve the likelihood that 
the aerosol sampled Is pure maritime in origin; because of 
the larger uncertainty on MAN AODs than those from the 
on-iaiid AERONET sites, the range of permitted Angstrom 
exponents is extended to -™0.1 < a < 1.2 . Unfortunately, 
many points remain in coastal or continental outflow regions, 
such that some contamination by a non-maritirne compo- 
nent is likely. To reduce this, it is further required that the 
measurement be at least 5° from land. This leaves 135 po- 
tential cases for comparison (only 104 include AOD at 440 
urn). 

Next the recommended aerosol size distribution param- 
eters and refractive index as outlined above (?\/.f = 0.157 
pm, rv,c = 2.58 pm, erf ~ 0.50, erf = 0.72, fine mode 
m. ~ 1.415 — 0.002i, coarse mode m = 1.363 ” 3 x 10“^i) 
are used with the relationships Cv,{ — 0.0056 and Cv,c ~ 
0.015 + 0.0036u;s to predict the MAN AOD and a. Statis- 
tics of the resulting comparison are shown in Table 1 1 , and 
a .scatter plot in Figure 17. Correlation coefficients are be- 
tween 0.43 and 0.54 for spectral AOD, which, although low, 
are significant, and reflect the low range of AODs fis well as 
the high variability of aerosol volume within a single nar- 
row wind speed bin previously observed in Figure 9. The 
SMAD is of similar magnitude to the uncertainty in MAN 
AOD. Because of all these reasons, a is poorly- reproduced 
overall. The majority of predicted AODs are within the 
MAN uncertainty. The wind-speed reURionship tends to 
slightly overestimate the AOD for low MAN AODs, and un- 
derestimate for high man AODs. This could imply either a 
stronger wind-AOD relationship over the open ocean than at 
the AERONET sites, or that some of the higher-AOD MAN 
observations do have a residual non-marine component. An 
alternative could be differences between the aerosol proper- 
ties for coastal and open-ocean regions, although this cannot 
be assessed as the MAN record does not permit retrievals of 
aerosol size distribution. Performance is similar if the other 
relationships betwfien wind speed and volume from Section 
3.2.1 are applied instead. The main conclusion from this 
is to confirm that the wind speed alone is not able to pre- 
dict the instantaneous aerosol burden well over the world’s 
oceans, although it can provide a typical background value. 


6. Lidar ratios 

Lidar provide a useful tool for examining the vertical 
structure of aerosols and clouds, which is not readily ac- 
cessible from radiometers to the same extent. To convert 
the backseat ter ing measured by a lidar from a particular 
altitude range into an extinction (which can then be inte- 
grated vertically to determine the optical depth), the lidar 
ratio S, which is a function of aerosol type and wavelength, 
is required. This can be calculated as 
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where P{0), the scattering phase function, is typically nor- 
malised to integrate to either 1 or 47 t, and evo is the aerosol 
single sccittering albedo; i.e. S is simply the ratio of total 
extinction to backscatter. Wavelength-dependence is omit- 
ted in the above for clarity. As such for a given vertically- 
integrated backscatter, the calculated AOD is directly pro- 
portional to 5. Depending on tlie characteristics of an in- 
dividual lidar, S is either calculated from measured extinc- 
tion and backscatter (for example, Pedros et al, 2010), or 
prescribed as a function of aerosol type. In the latter situ- 
ation the choice of an appropriate 3 is therefore important 
for the accurate calculation of aerosol loading. The Cloud- 
Aerosol Lidar with Orthogonal Polarisation (CALIOP) sen- 
sor has flown as part of the A-Train satellite constellation 
since 2006, and measures backseat tered radia,tion at 532 nm 
and 1064 nm (as well as depolarization at 532 nm); the li- 
dar ratios used in the current processing are given by Omar 
et a.L [2009]. For unpolluted marine aerosol, these are 20 
at 532 run and 45 at 1064 nm. Table 3 of Cattrall et al. 
[2005j provides a summary of some results from the litera- 
ture, with mid-visible S for marine aerosol between 24 and 
39. Their results also show little spectral dependence. Ack- 
er rriann [1998] performed calculations based on the OPAC 
database, giving S in the range 25-30 at 532 nm and 40-60 at 
1064 nni for typical relative humidities. Muller et al [2007] 
summarise a set of field campaigns, in which S at 532 nm 
was found to be from 23-29 for marine aerosol. Pedros et al. 
[2010] obtained median values of S at 532 nm of 31, 36, and 
37 for air masses originating from different oceanic regions, 
although there may have been a local contribution to these 
results. Results from Lanai for the AERONET analysis of 
Smirnov et al. [2003a] give 34.5 at 500 nm and 37 at 1020 
nm . 

Lidar ratios from distributions obtained in this study are 
presented in Table 12 for the commonly- used wavelengths of 
532 nm and 1064 nm, and are in the range 25-35. Because 
the bimodal lognormal fits do not reproduce perfectly the 
retrieved AERONET size distributions, two sets of calcula- 
tions are presented. All of these assume the refractive index 
m = 1.415 — 0,002i (fine mode) and m == 1.363 ~ 3 x 10“^' f 
(coarse mode). The first set uses the bimodal lognormal 
distribution parameters for the recommended aerosol model 
(Section 4), together with bimodal lognormal fits for the dis- 
tributions binned by wind speed (lower part of Table 4). The 
second uses the averaged size distributions directly, rather 
than lognormal fits to them. This makes the assumptions 
that the volume of particles outside the range of the bins 
is negligible (supported by Figure 3), and that the cutoff 
radius between fine and coarse mode (to determine which 
refractive index to use) is 0.4 /im, which is close to the in- 
flection point in Figure 3. If the cutoff radius is changed in 
the region 0.3 < r < 0.5 pm, S changes by approximately 
1 %, so sensitivity to this assumption is small. Additionally, 
as the recommended model was determined by the weighted 
average of lognormal fit parameters (rather than a fit to a 
weighted spectrum), results for Lanai are included in this 
second case (see Figure 4), as the site with the best sam- 
pling. This latter method is expected to yield more accu- 
rate lidar ratios, as it uses the retrieved size distributions 
directly. 

Use of exact distributions rather than lognormal fits re- 
sult in lidar ratios higher by approximately 10% at both 
wavelengths. This difference is due to the fact that the 
size distributions are not perfectly bimodal lognormal. The 
results are in good agreement with the ranges of the 
previously-cited studies; they are, however, significantly 
higher at 532 nrn (lower at 1064 nm) than the values used in 
CALIOP processing (20 at 632 nm and 45 at 1064 nm; Ornar 
et al, 2009). This is. an important result as applying a lidar 
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ratio of 30 as opposed to 20 for CALIOP at 532 nm would 
increase the unpolluted marine aerosol loading by 50%; this 
would explain the relative low bias of CALIOP marine AODs 
at 532 run as compared to other datasets over the ocean 
(Ktliyanpilakkil and Meskhidze, 2011, Oo and Holz. 2011). 
The difference at 1064 nrn is of similar magnitude but oppo~ 
site sense (i.e. CALIOP lidar ratio around 50% larger than 
the results here suggest). There is a slight decrease of li- 
dar ratio with wind speed; the change is typically 10% or 
less from the avera.ge value. This provides an estimate of 
the error which would arise from tlie assumption of a wind- 
speed-independent lidar ratio. 

If m — 1.37 - O.OOli is used (as in Smirnov et ai, 2003a), 
S increasrjs by approximately 10% at 532 nm and less at 
1064 nrn. The highligiits the sensitivity to the assumed re- 
fractive index, Additionally, the AERONET retrievals do 
not provide information on the vertical profile of the aerosol, 
provided a column- integrated amount. Therefore a strong 
vertical inhomogeneity in particle number or size may lead 
to errors in the retrieved size distributions, and influence the 
calculated S. Conversely, if there is significant vertical in- 
homogeneity in the aerosol size distribution or composition, 
the assumption of vertically-constant 5 will be inappropri- 
ate when trying to estimate total extinction from a lidar. 

7. Conclusions 
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Appendix A: Relation between number 
and volume size distributions 

For individual lognormal components, the relationships 
between the volume and number distribution parameters 
niciy be calculated using Equations 2 and 6 by first noting 
that (dropping subscripted i) 


When aerosol size distributions retrieved at 11 island 
AERONET sites spread throughout global oceans are fil- 
tered to extract data likely representative of unpolluted ma- 
rine aerosol, the resulting size distributions are similar, with 
the chief differences between sites being in the total fine and 
coarse mode volumes. An aerosol model with size distri- 
bution parameters and refractive index shown in Table 13 
was found to be able to reconstruct the AERONET AOD 
with accuracy of order 0.01-0,02, if only the fine and coarse 
mode volumes are taken as input. This accuracy is similar 
to that of the AERONET AOD measurements themselves, 
and holds at most sites and wavelengths between 340 nm 
and 1020 mn. These parameters are therefore suggested for 
use in aerosol remote sensing algorithms to represent unpol- 
luted marine aerosol. 

Size distributions were found to have a dependence on the 
wind speed, with higher winds leading to an increased coarse 
mode total volume and volume mean radius. As the major- 
ity of the data were for wind speeds between 4 and 8 ms“^ 
however, the global average coarse-mode radius can be used 
in most situations. The fine mode was coiriparatively unaf- 
fected. The AOD and Angstrom exponent also showed an 
approximately linear relationship with wind speed. How- 
ever, correlations were poor unless fits were performed to 
binned data, underlying the fact that v/ind speed alone is 
a poor predictor of the marine aerosol burden. When the 
relationship between wind speed and coarse mode aerosol 
volume was used to predict AOD observed on MAN cruises, 
the data were comparatively unbiased and had a scatter 
similar to the uncertainty on ship- measured AOD, although 
with a poor correlation, because of the small range of AOD. 
Relative liumidity was also invesitaged, although poor sam- 
pling and potential for influence of eitlier transported dust 
or cloud contamination for the lowest and highest humidities 
limit the strength of any conclusions which can be drawn. 
For the 88% of the data within 60 < rh < 80% there was 
little change in size distribution parameters. Similarly, SST 
and columnar water vapor were not found to have a strong 
impact on aerosol size distribution parameters, within the 
ranges sampled. 

Lidar ratios calculated from the results in this work give 
results approximately 50 % higher than those used for oper- 
ational CALIOP algorithm at 532 nm, but more consistent 
with those from many other studies, consistent with an un- 
derestimate of marine AOD by CALIOP as observed in other 
studies. 
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then expanding the exponential term 
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multiplying Equation A1 by 


_ g31n(fn) T 4.5cj 31n(rii) — 4.5(7^ (A3) 

and combining the second exponential factor in Equation 
A3 with Equation A2 to give 


“ [In^ (r) - 2 lii(rn ) ln(r) + In^ (rn) 

~~ 6(T^ in(r) f In(rn) 4- 9cr^j, (A4) 

which simplifies to 
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From this the definition 
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leads to the result that 
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Table 1, Locations of the AERONET sites used in this 
work, as well as the number and proportion of AERONET 
inversions passing the criteria for clean maritime conditions, 
as described in the text. 


Site name 

Latitude, 

degrees 

Longitude, 

degrees 

Elevation, m 
above sea level 

Number of 
maritime 

Proportion 
of maritime 

Lanai 

20.7350 

-156.922 

20 

1113 

0.79 

Bermuda 

32.3700 

-64.6960 

10 

116 

0.49 

Kaashidhoo 

4.96500 

73.4660 

0 

50 

0.20 

Midway Island 

28.2100 

-177.378 

20 

484 

0.86 

Ascension Island 

-7.97600 

-14.4150 

30 

341 

0.61 

Tahiti 

-17.5770 

-149.606 

98 

375 

0.82 

Amsterdam Island 

-37.8100 

77.5730 

30 

32 

0.68 

Crozet Island 

-46.4350 

51.8500 

221 

8 

0.47 

Guam 

13.4310 

144.801 

62 

74 

0.82 

Nauru 

-0.521000 

166.916 

7 

101 

0.91 

Graciosa 

39.0910 

-28.0300 

15 

265 

0.78 



X - 18 


SAVER ET AL.: PURE MARINE AEROSOL MODEL 


Table 2. Size distribution parameters for each site. The up- 
per half of the table shows averaged (median) size distril>ution 
parameter.s for AERO NET aerosol volume .size distributions, 
and figures in parentheses indicate the SMAD (defined in the 
text). The lower half shows bimodal lognormal distribution 
parameters for fits to averaged AERONET aeroscd volume size 
distributions, and figures in parentheses indicate one standard 
deviation uncertainty on the fit. Also shown for both cases are 
the mean values over all sites, and the mean weighted by the 
number of AERONET retrievals at each .site. 


Site name 

6\a, " 


■rv,f. fiin /Ain 


(7f 


Lanai 

0.0050 (0.0030) 

AER.ONET average parameters 
0.032 (0.013) 0.169 (0.019) 2.39 (0.25) 

0.4S 

(0.037) 

0.67 (0.040) 

Bermuda 

0.0080 (0.0044) 

0.041 (0.024) 

0.159 (0.024) 2.36 (0.51) 

0.46 

(0.0,36) 

0.65 (0,047) 

Kaashidhoo 

0.0080 (0.0030) 

0.037 (0.024) 

0.182 (0.013) 2.35 (0.47) 

0.45 

(0.043) 

0.71 (0.059) 

Midway Island 

0.0060 (0.0030) 

0.044 (0.024) 

0.167 (0.021) 2.41 (0.34) 

0.47 

(0.044) 

0.66 (0.042) 

Ascension Island 

0.0090 (0.0044) 

0.049 (0.019) 

0.156 (0.01,5) 2.36 (0.32) 

0.48 

(0.053) 

0.69 (0.033) 

Tahiti 

0.0040 (0.0015) 

0.028 (0.013) 

0.171 (0.021) 2.43 (0.33) 

0.48 

(0.034) 

0.69 (0.033) 

Amsterdam Island 

0.0050 (0.0030) 

0.028 (0.021) 

0.183 (0.0,34) 2.30 (0.,30j 

0.52 

(0.044) 

0.70 (0,0.39) 

Crozet Island 

0.0030 (0.0030) 

0.019 (0.022) 

0,251 (0.065) 2.17 (0.31) 

0.54 

(0.099) 

0,72 (0.077) 

Guam 

0.0060 (0.0030) 

0.035 (0.019) 

0.177 (0.031) 2.45 (0.29) 

0.49 

(0.031) 

0.69 (0.034) 

Nauru 

0,0040 (0.0015) 

0.031 (0.018) 

0.181 (0.031) 2., 55 (0.36) 

0.49 

(0.044) 

0.67 (0.044) 

Graciosa 

0.0060 (0.0030) 

0.030 (0.019) 

0.173 (0.018) 2..34 (0.39) 

0.52 

(0.049) 

0.74 (0.037) 

Mean 

0.0058 

0.034 

0.179 2.37 

0.49 

0.69 

Weighted mean 

0.0057 

0.036 

0.169 2.39 

0.48 

0.68 

Lanai 

0.0051 (0.0004) 

Bimodal fit to median distribution 

0.031 (0.002) 0.156 (0.007) 2.55 (0.11) 

0,50 

(0.049) 

0,72 (0.042) 

Bermuda 

0.0081 (0.0006) 

0.041 (0.002) 

0.145 (0.005) 2.54 (0.12) 

0.46 

(0.037) 

0.73 (0.049) 

Kaashidhoo 

0.0078 (0.0005) 

0.037 (0,002) 

0.170 (0.005) 2.56 (0.13) 

0.44 

(0.032) 

0.76 (0.052) 

Midway Island 

0.0056 (0.0004) 

0.043 (0.002) 

0.157 (0.006) 2.58 (0.10) 

0.49 

(0.043) 

0.70 (0.039) 

Ascension Island 

0.0082 (0.0008) 

0,047 (0.003) 

0.139 (0.007) 2.,56 (0.1.3) 

0.45 

(0.054) 

0,73 (0.050) 

Tahiti 

0.0(,)41 (0.0003) 

0,027 (0.002) 

0.161 (0.007) 2.72 (0.14) 

0.51 

(0.047) 

0.69 (0.0,50) 

Amsterdam Island 

0.0037 (0.0002) 

0.030 (0.002) 

0.155 (0.004) 2.29 (0.15) 

0.61 

(0.030) 

0.84 (0.066) 

Crozet Island 

0.00.30 (0.0002) 

0,022 (0.002) 

0.190 (0.006) 2.54 (0.25) 

0.49 

(0.034) 

0.94 (0.10) 

Guam 

0.0054 (0.0004) 

0.034 (0.002) 

0,170 (0.008) 2.62 (0.10) 

0.56 

(0.053) 

0.74 (0.039) 

Nauru 

0.0035 (0.0003) 

0,031 (0.002) 

0.174 (0.009) 2.94 (0.14) 

0.55 

(0.056) 

0.67 (0.049) 

Graciosa 

0.0062 (0.0005) 

0,029 (0,002) 

0,168 (0.008) 2.48 (0.16) 

0.57 

(0.052) 

0.83 (0.068) 

Mean 

0.0055 

0.034 

0.162 2.58 

0.50 

0.76 

Weighted Mean 

0.0056 

0.035 

0.157 2.58 

0.50 

0.72 
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Table 3. Number of retrievals, AOD at two wavelengths, 
and Angstrom exponent corresponding to AERO NET aerosol 
volume size distributions, binned as a function of near-surface 
wind speed. Figures in parentheses indicate the SMAD. 


Wind 

Number of 

TA-%0 

no2o 


a 

speed 

retrievals 







0-4 

437 

0.068 

(0.022) 

0.040 

(0.013) 

0.73 

(0.14) 

4-6 

1022 

0.065 

(0.019) 

0.039 

(0.013) 

0.60 

(0.17) 

6-8 

1082 

0.076 

(0.025) 

0.052 

(0.018) 

0.49 

(0.18) 

8-10 ms~^ 

299 

0.081 

(0.025) 

0.058 

(0.020) 

0.37 

(0.15) 

10+ ms~^ 

67 

0.090 

(0.037) 

0.062 

(0.02.1) 

0.34 

(0.14) 
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Table 4. Size distribution parameters for average distribu- 
tions binned as a function of near-surface wind speed. The 
upper half of the table shows a,veraged (median) size distribu- 
tion parameters for AERONET aerosol volume size distribu- 
tions, and figures in parentheses indicate the SMAD (defined 
in the text). The lower half shows bimodal lognormal dis- 
tribution parameters for fits to averaged AERONE'F aerosol 
volume size distributions, and figures in parentheses indicate 
one standard deviation uncertainty on the fit. 


Wind spec'd 


—p; T! 

Ov,c fim ftm 

rv,f, fim 

rv,c. /ira 


<Tc 

refrj, /.nn 

reff,c Rm 




AERONET average parameters 




0-4 ms“ ^ 

0.0060 (0.0020) 

0,024 (0.0090) 

0.173 (0.013) 

2.27 (0.25) 

0.46 (0.027) 

0.70 (0.041) 

0.155 (0.012) 

1.76 (0.18) 

4-6 ms ^ 

0,0060 (0.0020) 

0.030 (0.0090) 

0.168 b-014) 

2.38 (0.20) 

0.48 (0.026) 

0,69 (0-029) 

0.150 (0.012) 

1.87 (0.14) 

6-8 ms '■ 

0.0060 b-0020) 

0,040 (0.012) 

0.167 b-014) 

2.43 b-19) 

0.49 (0.026) 

0.68 (0.027) 

0.148 5-012) 

1.90 (0.14) 

8-10 

0.0060 (0.0020) 

0.051 (0-016) 

0.168 b-014) 

2.45 (0-19) 

0.50 (0.028) 

0.68 b-031) 

0.149 (0.011) 

1.92 (0.15) 

10+ ms“^ 

0.00.50 (0.0020) 

0.054 b-021) 

0,169 (0.020) 

2.45 (0-28) 

0.52 (0,018) 

0.70 (0.040) 

0.147 (0-016) 

1.91 5-18) 




Bimodal fit to i 

median distribution 




0-4 ms"^ 

0.0059 (0,00034) 

0.024 (0.0012) 

0.167 (0.0054) 

2.35 (0.11) 

0.50 (0.035) 

0.79 (0.047) 

0.131 (0.0042) 

1.58 (0.073) 

4-6 ms 

0.0052 (0.00040) 

0.030 (0.0016) 

0.156 (0.0067) 

2.56 (0.12) 

0.49 (0.046) 

0.73 (0.046) 

0.122 (0-0052) 

1.77 (0.080) 

6~8 

0.0055 (0.(»050) 

0.040 (0.0022) 

0.152 (0.0079) 

2.63 (0.12) 

0.51 (0.056) 

0.71 (0.046) 

0.118 (0.0061) 

1.84 (0.084) 

8-10 ms ’ ^ 

0.0055 b-00053) 

0.048 (0.0028) 

0.155 b-0089) 

2.72 b-13) 

0.53 (0.062) 

0.69 (0.048) 

0.119 5-0068) 

1.92 5-091) 

10+ ms"^ 

0.0044 (0.00033) 

0.0-49 03.0030) 

0.143 ii.ooeo) 

2.70 (0.14) 

0.50 b-048) 

0.76 (0.054) 

0.111 il.0047) 

1.84 5-099) 


o 



Table 5. Statistics of linear regression between wind speed 
and AOD or a (daily averages for both datasets) of the form 
r (or q)~ a + 6 X ws, from direct-Sun and SDA AERONET 
data. Also shown is Pearson’s linear correlation coefficient for 
the fit, R. The bottom row shows the fits when data from all 
sites are considered together. 


Site 

a 

T440 

b 

i? 

a 

T 1020 

h 

R 

a 

Tf .500 

b 

R 

tt 

Tc, 500 

h 

R 

G: 

a b 

R 

Lanai 

0.092 

-0.0021 

0.13 

0.044 

0.00058 

0.053 

0.082 

-0.0015 

0.10 

0,042 

-0.0021 

0.23 

0.79 -0.039 

0.35 

Bermuda. 

0.097 

0.0028 

0.20 

0.053 

0.0029 

0.29 

0.088 

0.0026 

0.20 

0.046 

-0.000021 

0.0031 

0.82 -0.024 

0.31 

Kaa»shidoo 

0.12 

0.0022 

0.091 

0-070 

0.0033 

0.18 

0.11 

0.0023 

0.10 

0.052 

-0.0013 

0.11 

0.70 -0-036 

0.27 

Midway Island 

0.069 

0.0044 

0.25 

0.033 

0.0051 

0.37 

0.063 

0.0044 

0.27 

0.034 

-0.00017 

0.023 

0.76 -0.039 

0.40 

Ascension. Island 

0.090 

0.0053 

0.13 

0.047 

0.0039 

0.12 

0-080 

0.0052 

0.14 

0-039 

0.00070 

0.034 

0.71 -0.024 

0.11 

Taiiiti 

0.082 

0.00041 

0,024 

0.029 

0.001 1 

0.058 

0.075 

0.00043 

0.026 

0.032 

-0.00014 

0.020 

0.72 -0.021 

0.17 

Amsterdam Island 

0.057 

0.0030 

0.30 

0.033 

0.0032 

0.37 

0.051 

0.0030 

0.32 

0.016 

0.0011 

0.27 

0.58 -0.013 

0.17 

Cr 02 ,et Island 

0-061 

0.0015 

0.12 

0.044 

0.0016 

0.16 

0.054 

0.0017 

0.14 

0.021 

0-00052 

0.072 

0.55 -0.011 

0.13 

Guam 

0.097 

0.000070 

0.0029 

0.051 

0.0031 

0.16 

0.087 

0.0011 

0.048 

0.051 

-0.0023 

0.17 

0.90 -0.074 

0.48 

Naum 

0.055 

0.0041 

0.24 

0.028 

0.0047 

0.23 

0.052 

0.0043 

0.26 

0.022 

0.00055 

0.087 

0.60 -0.039 

0.26 

Graciosa 

0.084 

0.0051 

0.31 

0.047 

0.0061 

0-45 

0.077 

0.0053 

0-34 

0.044 

-0-00049 

0.074 

0.75 -0.044 

0.45 

All sites 

0.071 

0.0036 

0.20 

0-033 

0.0040 

0.27 

0-066 

0.0034 

0.21 

0.031 

0.00038 

0.048 

0.65 -0.022 

0.21 


X 

to 
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Table 6. Number of retrievals, AOD at two wavelengths, 
Angstrbrn exp{:>neiit, and average wind speed corresponding 
to AERO NET aerosol volume size distributions, binned as a 
function of near-surface relative humidity. Figures in paren- 
theses indicate the SMAD. 


Relative 

humidity 

jM urn her of 
retrievals 

T440 

T1020 

a 

Wind 

speed 

0 - 50 % 

4 

0.10 (0.021) 

0.082 (0.017) 

0.22 (0.10) 

5.07 

50 - 60% 

67 

0.081 (0.027) 

0.058 (0.014) 

0.35 (0.18) 

6.42 ms” ' 

60 ■ 65 % 

270 

0.068 (0.020) 

0.044 (0.014) 

0.49 (0.21) 

6.53 

65 - 70% 

718 

0.068 (0.02.3) 

0.044 (0,016) 

0.49 (0.17) 

6.77 ms-‘ 

70 - 75 % 

920 

0.074 (0.024) 

0.045 (0.017) 

0.59 (0.18) 

6.13 rns“' 

75 - 80% 

647 

0.074 (0.022) 

0.046 (0.016) 

0.60 (0.19) 

5.53 ms'-‘ 

80 - 85 % 

237 

0,069 (0.019) 

0.044 (0.014) 

0.66 (0.15) 

5.10 m.s~* 

85 - 90% 

38 

0.060 (0.011) 

0.050 (0.0098) 

0.41 (0.14) 

5.56 m.s-' 

90 - 100% 

4 

0.16 (0.0052) 

0,14 (0.0068) 

0.21 (0.021) 

6.44 ms"‘ 



lable 7. As Table 4, except for AERONET distributions binned as a function of relative humidity. 

Relative Cv,f, Cvx, 7’v,f, /im rv,c, Mui ijf 

humidity 


0-50% 0.0050 (0.0010) 0.073 (0.0090) 

50 - 60% 0.0050 (0.0020) 0.051 (0.012) 

60 - 65% 0.0050 (0.0020) 0.035 (0.012) 

65 - 70 % 0.0050 (0.0020) 0.036 (0.011) 

70 - 75 % 0.0060 (0.0020) 0.035 (0.012) 

75 - 80% 0.0060 (0.0020) 0.033 (0.012) 

SO - 85 % 0.0060 (0.0020) 0.031 (0.013) 

85 - 90% 0.0050 (0.0020) 0.032 (0.011) 

90- 100 % 0.0080 (0.0010) 0.13 (0.023) 

0-50% 0.0049 (0.00037) 0.070 (0.0050) 

50 - 60% 0.0050 (0.00047) 0.048 (0.0029) 

60 - 65% 0.0051 (0.00039) 0.035 (0.0019) 

65 - 70% 0.0049 (0.00046) 0.036 (0.0019) 

70 - 75 % 0.0057 (0.00047) 0.034 (0.0019) 

75 - 80 % 0.0058 (0.00044) 0.032 (0.0017) 

80 - 85% 0.0054 (0,00030) 0.030 (0.0017) 

85 - 90% 0,0043 (0.00042) 0.032 (0.0021) 

90 - 100 % 0.0087 (0.00044) 0.13 (0.0099) 


AERONET average parameters 
0,186 (0.026) 2.49 (0,088) 0.53 (0.013) 

0.169 (0.017) 2.51 (0.20) 0.49 (0.031) 

0.173 (0.014) 2.40 (0.20) 0.49 (0.027) 

0,168 (0.014) 2.46 (0.18) 0.49 (0.026) 

0.166 (0.014) 2.38 (0.19) 0.48 (0.026) 

0.168 (0.014) 2.33 (0.23) 0.48 (0.029) 

0.171 (0.013) 2.37 (0.28) 0.47 (0.027) 

0.174 (0.013) 2.21 (0.20) 0.54 (0.043) 

0.203 (0.032) 2.57 (0.24) 0.53 (0.00010) 

Bimodal fit to median distribution 
0.156 (0.0067) 2.89 (0.15) 0.51 (0.049) 

0.159 (0.0091) 2.86 (0.14) 0.56 (0.064) 

0.156 (0.015) 2.58 (0.12) 0.50 (0.027) 

0.153 (0.0083) 2.66 (0.11) 0.52 (0.059) 

0.154 (0.0072) 2.59 (0.12) 0.50 (0.050) 

0.156 (0.0065) 2.47 (0.11) 0.49 (0.045) 

0.161 (0.0050) 2.55 (0.13) 0.49 (0.033) 

0.155 (0.0093) 2.33 (0.14) 0.56 (0,067) 

0.214 (0.0081) 2.68 (0.20) 0.67 (0.041) 


«rc reff.f, lira Tch.c, /rm 


0.64 (0.025) 0.162 (0.022) 1.99 (0.056) 

0.68 (0.027) 0.150 (0.014) 1.98 (0.15) 

0.67 (0.028) 0.154 (0.011) 1.90 (0.15) 

0.68 (0.027) 0.149 (0.013) 1.93 (0.14) 

0.69 (0.028) 0.149 (0.012) 1.87 (0.14) 

0.69 (0.035) 0.150 (0.012) 1.82 (0.15) 

0.69 (0.036) 0.152 (0.012) 1.84 (0.19) 

0.72 (0.018) 0.152 (0.0090) 1.71 (0.15) 

0.74 (0.012) 0.175 (0.027) 1.96 (0.19) 

0.62 (0.051) 0.121 (0.0052) 2.13 (0.11) 

0.68 (0.048) 0.121 (0.0069) 2,03 (0.097) 
0.72 (0.045) 0.122 (0.012) 1.80 (0.081) 
0.71 (0.043) 0.118 (0.0064) 1,86 (0.080) 
0.72 (0.046) 0.120 (0.0056) 1.80 (0.082) 
0.75 (0.045) 0.122 (0.0051) 1.70 (0.075) 
0.77 (0.051) 0.126 (0.0039) 1.74 (0.087) 

0-83 (0.064) 0.117 (0.0070) 1.53 (0.094) 
0.87 (0.080) 0.153 (0.0058) 1.74 (0.13) 
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Table 8 . Number of retrievals, AOD at two wavelengths. 
Angstrom exponent, and average wind speed corresponding 
to AERONET aerosol volume size distributions, binned as 
a function of AERONET columnar water vapor. Figures in 
parentheses indicate the SMAD. 

Water Number of T 440 ^1020 Wind 

vapor retrievals speed 


0~1 cm 

15 

0.034 

(0.019) 

0.029 

(0.020) 

0.35 

(0.15) 

5.30 

ms' ^ 

L2 cm 

508 

0.077 

(0.027) 

0.052 

(0.019) 

0.46 

(0.20) 

6.44 

ms-‘ 

2-2.5 cm 

615 

0.065 

(0.021) 

0.043 

(0.01.5) 

0.51 

(0.18) 

6.56 

ms'”^ 

2.5-3 crn 

713 

0.070 

(0.022) 

0.044 

(0.015) 

0.60 

(0.19) 

6.00 

ms""^ 

3-3.5 cm 

620 

0.074 

(0.021) 

0.046 

(0.016) 

0.60 

(0.18) 

5.82 


3.5-4 cm 

291 

0.071 

(0.018) 

0.046 

(0.016) 

0.62 

(0.16) 

5.91 


4+ crn 

143 

0.071 

(0.019) 

0.043 

(0.014) 

0.63 

(0.20) 

5.27 
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Table 9. Refractive indices m ~ n~~ ik for marine aerosol tested in this work, along with information on sources. 


Ciise 

Fine mode 

Coarse mode 

Note 

1 

1.37 - O.OOli 

1.37 - O.OOli 

Smirnov et al [2003al; average of AERONET 
inversion results for Lanai, 

2 

1.45 - 0.0035i 

1.35 - O.OOli 

Rerner et al. [2006]; pair of oceanic components 
used in MODIS collection 5 aerosol retrieval. 

3 

1.39 - 0.003i 

1.39 » 0.0031 

Silva et ai [2002]; ground-based unpolluted 
maritime measurements on the Portuguese coast. 

4 

L415 - 0.G02i 

1.363 - 3 X 10 

Shettle and Fenn, 1979/ //ess et al, 1998 at 
500 nm and rh —70 %; fine: water soluble 
component; coarse: accumulation /coarse sea 
s alt component . 

5 

1.415 - 0.002i 

1.434 - 3 X 10”^i 

As CJise 4j except real part of coarse mode from 
Irshad et al [2009] infrared laboratory data 
at rh —74.2%, retrieval value at 500 nm. 
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Table 10. Pearson's linear correlation coefficient (top por- 
tiori), .median (calculated ~ AERONET observed) bias (mid- 
dle portion), and SMAD (bottom portion) between observed 
AERONET spectral AOD and Angstrom exponent, and that 
calculated using the average ‘lognormal fitted’ fine and coarse- 
mode radii and spreads, for refractive index case 4 from Table 
9. 


Site 

r;,M0 

T:m 

T440 

"500 

T675 

r B70 

TI.020 

a 

Lanai 

0.95 

Correlation coefficient 
0.96 0.97 0.96 

0.95 

0.94 

0.94 

0.77 

Bermuda 

0.91 

0.91 

0.92 

0.91 

0.88 

0,87 

0.87 

0.77 

Kaashidhoo 

0.89 

0.90 

0,91 

0.90 

0.89 

0.89 

0.88 

0.88 

Midway Island 

0.90 

0.91 

0.92 

0.92 

0.92 

0.93 

0.94 

0.88 

Ascension Island 

0.88 

0.89 

0.89 

0.89 

0.88 

0.88 

0.89 

0.88 

Tahiti 

0.89 

0.91 

0.95 

0.94 

0,94 

0.95 

0.93 

0.63 

Amsterdam Island 

0.95 

0.95 

0.96 

0.96 

0.96 

0.96 

0.97 

0.74 

Crozet Island 

0.96 

0.95 

0.93 

0.93 

0.85 

0.67 

0.68 

-0.97 

Guam 

0.84 

0.84 

0.84 

0.84 

0.83 

0.79 

0.78 

0.76 

Nauru 

0.78 

0.82 

0.85 

0.86 

0.91 

0,90 

0.90 

0.79 

Graciosa 

0.76 

0.80 

0.82 

0.84 

0.89 

0.90 

0.91 

0.72 

Lanai 

0.0015 

0.0014 

Median bias 
0.00067 -0.0026 

-0.0015 

-0.0052 

-0,00030 

0.17 

Bermuda 

-0.0094 

-0.011 

-0,0094 

-0.010 

-0.0086 

-0.0090 

-0.0095 

0.12 

Kaashidhoo 

-0.033 

-0.031 

-0.023 

-0.025 

-0.022 

-0.018 

-0.019 

0.15 

Midway Island 

-0.00045 

-0,0048 

-0.0011 

-0.0045 

-0.0068 

-0.0086 

-0.0015 

0.22 

Ascension Island 

-0.0077 

-0.0090 

-0.0071 

-0.011 

-0.011 

-0.014 

-0.0095 

0.19 

Tahiti 

-0.012 

-0.0042 

-0.0040 

-0.0065 

-0.0029 

-0.0041 

-0.0029 

0.050 

Amsterdam Island 

0.0068 

-0.00042 

-0.00097 

-0.0037 

-0.0042 

-0.0099 

-0.0068 

0.20 

Crozet Island 

-0.0068 

-0.0063 

-0.0091 

-0.0068 

-0.0041 

-0.012 

-0.014 

0.74 

Guam 

0.0064 

0.0047 

0.0029 

-0,0029 

-0.0016 

-0.0056 

-0.0026 

0.15 

Nauru 

-0.0015 

-0.0033 

-0.0017 

-0.0046 

-0.0041 

-0.0042 

-0.0043 

0.096 

Graciosa 

-0.021 

-0.013 

-0.012 

-0.018 

-0.013 

-0.013 

-0.016 

0.18 

Lanai 

0.011 

Scaled median absolute difference 
0.0081 0.0057 0.0058 0.0061 

0.0056 

0.0061 

0.13 

Bermuda 

0.021 

0.021 

0.018 

0.018 

0.016 

0.015 

0.014 

0.19 

Kaashidhoo 

0.013 

0.0094 

0.0062 

0.011 

0.0068 

0.0060 

0.0056 

0.19 

Midway Island 

0.014 

0.013 

0.010 

0.0099 

0.0084 

0.0078 

0.0078 

0.10 

Ascension Island 

0.018 

0.016 

0.016 

0.016 

0.015 

0.012 

0.011 

0.11 

Tahiti 

0.011 

0.0074 

0.0050 

0.0053 

0-0049 

0.0046 

0.0059 

0.15 

Amsterdam Island 

0.011 

0.0094 

0.0080 

0.0081 

0.0091 

0.0082 

0.0064 

0.10 

Crozet Island 

0.025 

0.025 

0.021 

0.028 

0.027 

0,012 

0.022 

0.80 

Guam 

0.010 

0.011 

0.0087 

0.0072 

0.0059 

0,0065 

0.0079 

0.12 

Niiuru 

0.0086 

0.0089 

0.0071 

0.0089 

0.0049 

0.0076 

0.0074 

0.13 

Graciosa 

0.017 

0.014 

0,013 

0.011 

0.0093 

0,0093 

0.0087 

0.16 
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Table 11. Comparison between spectral AOD and a mod- 
elled using NCEP wiiid speed, and that measured on MAN 
cruises. The correlation coefficient i.s Pearson ts linear corre- 
lation coefficient, and the bias is defined such that positive 
values mean the model is larger than the MAN data. The fi- 
nal column shows, for spectral AOD, the fraction of predicted 
AODs lying within the Sun- photometer typical uncertainty of 
± 0 . 02 . 


Parameter 

Number of 
matches 

Correlation 

coefficient 

Median 

bias 

SMAD 

AHnimuin 
MAN value 

Maximum 
MAN value 

Fraction 
within 0.02 

T440 

104 

0.45 

0.001 

0.024 

0.022 

0.235 

0.55 

TyOD 

135 

0.43 

0.006 

0.021 

0,021 

0.175 

0.62 

T675 

135 

0.46 

0.005 

0,020 

0.017 

0.143 

0.68 

T870 

135 

0.54 

0.0002 

0.018 

0.017 

0.143 

0.68 

a 

104 

0.23 

0.25 

0.322 

-0.053 

1.16 

- 
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Table 12. Lidar ratios S for unpolluted marine aerosol, cal- 
culated for 532 rim and 1064 nm. The top section presents 
results for bimodal lognormal distribution fits. The bottom 
section presents calculations for the exact siveraged size distri- 
butions. 


Distribution 

Lidar n 
532 mix 

atio S 
1064 nm 

Bimodal lognormal fits 


R,ecom mended model 

28.1 

30.8 

0-1 ms"' 

34.0 

35.2 

4-6 ms"^ 

29.3 

31.2 

6-8 rns”^ 

26.9 

29.9 

8-10 ms'^^ 

25.4 

28.7 

10+ ms" ^ 

25.0 

30.7 

Exact size distributions 


Lanai 

30.3 

33.3 

0-4 ms'"* 

35.9 

35.9 

4-6 ms“^ 

31.3 

33.6 

6-8 ms“^ 

29.3 

33.0 

8-10 ms"^ 

28.2 

32.3 

10+ms“^ 

28.4 

34.2 
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Table 13. Parameters for recommended aerosol model of 
unpolluted marine aerosoL for use in general satellite remote 
sensing applications. 

Parameter Value 
Fine mode 
7\rj 0,157 

o-r 0.50 

rn 1.415 - 0.002t 

Coarse mode 
?%-,c 2.58 fim 

(7c 0.72 

m 1.363 - 3 X 
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Locations of AERONET sites used 



Longitude, degrees 


Figure 1. Locations of AERONET sites used in this work. 
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Midway Island 


Kaashidhoo 




Figure 2. Averaged (median) retrieved aerosol volume 
size distributions at Midway Island (left) and Kaashidhoo 
(right), calculated from all retrieved size distributions 
(black), and only those distributions meeting the mar- 
itime criteria described in the text (red). The shaded re- 
gion indicates the to 95*^^ percentiles of dV{r)/dln{r) 
encountered in each size bin. 
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Figure 3. Averaged (median) aerosol volume size distri- 
butions for the AERONET sites in Table 1, for measure- 
ments corresponding to maritime- type size distributions 
as described in the text. 
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0.02 


O.Ol 
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Figure 4. Averaged (inedian) retrieved aerosol volume 
size distributions for Lanai (black), and lognormal ap- 
proximations to it. The distribution constructed from 
direct use of AERO NET parameters is in red, the best- 
fit bimodal lognormal distribution is in green, and the 
best-fit trimodal lognormal distribution in blue. Error 
bars on the retrieved size distribution denote the scaled 
median absolute deviation, as described in the text. 
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InstanLincous data 



0 5 10 J5 20 25 30_ 35 

Shif>base,d wind spt:ed, ms'" 


Daily average data 
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Ship-based wind speed, ms ' 


Instantaneous data 



20 40 60 80 100 

Ship-based relative humidity, % 


Figure 5, Scatter-density comparison between NCEP 
and ship- based measurements of wind speed and relative 
humidity. From top-bottom, plots show instantaneous 
wind speeds; daily- aver aged wind- speeds; and instanta- 
neous relative humidity. Bins without data are shown in 
white. The bin size is 0.5 ras" '^ for wind speed and 2% 
for relative humidity. The 1:1 line is overplotted. 
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Figure 6. Pearson’s linear correlation coefficient be- 
tween time- averaged wind speed and AERO NET aerosol 
properties, for all data combined (black lines) and those 
sites with 100 or more AERONET retrievals (coloured 
lines). Data for the AERONET aerosol size di.stributioii 
product (Dubovik and King, 2000). 
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0.04 
0.03 
0.02 
o.m 
0.00 

Figure 7. Averaged (median) AERONET aerosol volume size distributions, binned by near-surface wind speed. 
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Figure 8. Linear (red) and exponential (green) re- 
gression fits of aerosol volume to 24-hour- averaged wind 
speed for fine (left) and coarse (right) mode data from all 
AERONET sites considered together. Regression statis- 
tics are given in the plots, R is Pearson’s linear correla- 
tion coefficient for the fit. 
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Figure 9. Relationship between wind speed and aerosol 
volume for fine (left) and coarse (right) modes, binned 
by wind speed in bins of 0.5 ms~L Black diamonds show 
data binned by spatio-teinporally interpolated NCEP 
wind speeds, and red triangles data binned by NCEP 
wind speed averaged over the 24~hour~period prior to the 
retrieval. Error bars show the standard deviation on each 
bin's data. Coefficients of linear fit are given in the plots, 
and illustrated wdth dashed lines; R indicates Pearson’s 
linear correlation coefficient. Data are only shown where 
a bin contains at least 10 data points. 
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Figure 10. Relationship between daily averages of wind 
speed and AOD at 440 nrn (left), AOD at 1020 nra (mid- 
dle), and a (right), binned by wind speed in bins of 0.5 
ms”D Bin medians are used, and error bars show the 
standard deviation on each bin's data. Coefficients of lin- 
ear fit are given in the plots, and illustrated with dashed 
lines; R indicates Pearson’s linear correlation coefficient. 
Data are only shown where a bin contains at least 10 data 
points. 
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Figure 11. Averaged (median) AEROTsfET aerosol 
volume size distributions, arranged by near-surface 24- 
hour-averaged wind speed (top to bottom, range indi- 
cated above plots) , and binned according to the 24-hour- 
averaged SST, The number of size distribution inversions 
contributing to each line is given in the legend. 
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Figure 12. As Figure 7, except binned by relative humidity. 
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Figure 13. Relationship between daily averages of rel- 
ative humidity and AOD at 440 nm (left), AOD at 1020 
niri (middle), and a (right), binned by relative humidity 
in bins of 2 %. Bin medians are used, and error bars show 
the standard deviation on each bin’s data. Coefficients 
of linear fit are given in the plots, and illustrated v/ith 
dashed lines; R indicates Pearson’s linear correlation co- 
efficient. Data are only shown where a bin contains at 
least 10 data points. 




dV/dlnr, |inr|um'" 


SAYER ET AL,: PURE AIARINE AEROSOL MODEL 


X- 43 


0.03 


0.02 F 


0.0 1 F 


0.00 


Figure 14. 


Size distributions by water vapor 



As Figure 7, except binned by AERONET columnar water vapor amount. 
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Figure 15. Relationship between water vapor con- 
tent (wv) and r /,500 (left) or Tc,mo (right) from the 
AERO NET SDA product, for four sites. In each case, the 
red line indicates an exponential fit to all points, and the 
green line an exponential fit to only those points where 
T500 < 0.2 and 0.1 < a < 1 (deemed maritime condi- 
tions). The equation of each fit line, and Pearson’s linear 
correlation coelRcient of the fit (/?), are given in each 
plot. 
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Figure 16. Statistics of comparison of spectral AOD 
and Angstrom exponent between AERONET retrievals, 
and calculations performed using the average aerosol size 
distribution parameters, for a variety of assumed aerosol 
refractive indices (cases in Table 9). Subfigures show 
Pearson’s linear correlation coefficient (left), median bias 
(centre), and SMAD (right). Red diamonds indicate the 
‘AERONET lognormal ’ approach and green triangles the 
‘lognormal fitted’ method, slightly offset along the ;r-axis 
for clarity. Symbols show the medifin, and error bars the 
minimum and maximum values, over the ensemble of five 
sites used for the evaluation. 
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Figure 17. Comparison between MAN AODs and 
those predicted using the wind-speed relationship from 
AERONET sites. The solid line is the 1:1 line, and the 
dotted lilies indicate the MAN uncertainty of ±01)2. 




